Proton- Antiproton Annihilation 
and Meson Spectroscopy 
with the Crystal Barrel* 



OS 

os 
< 

00 



> 

in 
o 

oo 
O 

i> 

Os 

">< 
<D 

Oh! 



X 



Claude Amsler 
Physik-Institut der Universitat Zurich, 
Winterthurerstrasse 190, CH-8057 Zurich, Switzerland 

February 5, 2008 



Abstract 

This report reviews the achievements of the Crystal Barrel experiment at the Low En- 
ergy Antiproton Ring (LEAR) at CERN. During seven years of operation Crystal Barrel 
has collected very large statistical samples in pp annihilation, especially at rest and with 
emphasis on final states with high neutral multiplicity. The measured rates for annihila- 
tion into various two-body channels and for electromagnetic processes have been used to 
test simple models for the annihilation mechanism based on the quark internal structure 
of hadrons. From three-body annihilations three scalar mesons, ao(1450), /o(1370) and 
/o(1500) have been established in various decay modes. One of them, /o(1500), may be 
identified with the expected ground state scalar glueball. 
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1 Introduction 

Low energy antiproton-proton annihilation at rest is a valuable tool to investigate phenom- 
ena in the low energy regime of Quantum Chromodynamics (QCD). Due to the absence 
of Pauli blocking, the antiproton and proton overlap and one expects the interactions 
between constituent quarks and antiquarks (annihilation, pair creation or rearrangement) 
to play an important role in the annihilation process. From bubble chamber experiments 
performed in the sixties (Armenteros and French, 1969) one knows that annihilation pro- 
ceeds through qq intermediate meson resonances. The w(782), /i(1285), £'/r/(1440) and 
-KTi(1270) mesons were discovered and numerous properties of other mesons (ao(980), 
K*(892), 0(1020), a2(1320)) were studied in low energy pp annihilation. With the advent 
of QCD one now also predicts states made exclusively of gluons (glueballs), of a mixture 
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of quarks and gluons (hybrids) and multiquark states, all of which can be produced in pp 
annihilation. 

With the invention of stochastic cooling and the operation of the Low Energy An- 
tiproton Ring (LEAR) from 1983 to 1996, intense and pure accelerator beams of low 
momentum antiprotons between 60 and 1940 MeV/s were available at CERN. It is im- 
pressive to compare the high flux of today's antiproton beams (> 10 6 p/s) with the rate 
of about 1 p every 15 minutes in the early work when the antiproton was discovered, back 
in 1955 (Chamberlain, 1955). 

This survey covers the results obtained with the Crystal Barrel, designed to study low 
energy pp annihilation with very high statistics, in particular annihilation into n charged 
particles (n-prong) and m neutrals (ir°, rj, rj or uj) with m > 2, leading to final states with 
several photons. These annihilation channels occur with a probability of about 50% and 
have not been investigated previously. They are often simpler to analyze due to C-parity 
conservation which limits the range of possible quantum numbers for the intermediate 
resonances and the pp initial states. 

The experiment started data taking in late 1989 and was completed in autumn 1996 
with the closure of LEAR. Most of the data analyzed so far were taken with stopping 
antiprotons in liquid hydrogen on which I shall therefore concentrate. This article is 
organized as follows: After a brief reminder of the physical processes involved when an- 
tiprotons are stopped in liquid hydrogen (section 2), I shall describe in section 3 the 
Crystal Barrel apparatus and its performances. The review then covers results relevant 
to the annihilation mechanism and the roles of quarks in the annihilation process (section 
Q). Electromagnetic processes are covered in section |5| The observation of a strangeness 
enhancement may possibly be related to the presence of strange quarks in the nucleon 
(section |^). After describing the mathematical tools for extracting masses and spins of 
intermediate resonances (section |^) I shall review in sections 8 to 10 what is considered 
to be the main achievement, the discovery of several new mesons, in particular a scalar 
(J p = 0+) state around 1500 MeV, which is generally interpreted as the ground state 
glueball. Section 11 finally describes the status of pseudoscalars in the 1400 MeV region. 

In this review I shall concentrate on results published by the Crystal Barrel Collabo- 
ration or submitted for publication before summer 1997[]. Alternative analyses have been 
performed by other groups using more flexible parametrizations and also data from pre- 
vious experiments (e.g. Bugg (1994, 1996), Abele (1996a)). I shall only refer to them 
without describing them in detail since they basically lead to the same results. Results on 
pd annihilation will not be reviewed here. They include the observation of the channels 
pd — > 7r°n, 7]n, urn (Amsler, 1995a), pd A(1232)vr° (Amsler, 1995b) which involve both 
nucleons in the reaction process. 

2 Proton- Antiproton Annihilation at Rest 

Earlier investigations of low energy pp annihilation have dealt mainly with final states in- 
volving charged mesons (tt^, K^) or K$ — ► tt + 7t~, with at most one missing (undetected) 
7T°, due to the lack of a good 7 detection facility (for reviews, see Armenteros and French 
(1969), Sedlak and Simak (1988) and Amsler and Myhrer (1991)). 

The average charged pion multiplicity is 3.0 ± 0.2 for annihilation at rest and the 
average ir° multiplicity is 2.0 ± 0.2. The fraction of purely neutral annihilations (mainly 
from channels like 37r°, 5tt°, 2tt°t] and 4tt°t] decaying to photons only) is (3.9 ± 0.3) 
% (Amsler, 1993a). This number is in good agreement with an earlier estimate from 

All Crystal Barrel publications are listed with their full titles in the reference section. 
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bubble chambers, (4.1 ^gg) % (Ghesquiere, 1974). In addition to pions, rj mesons are 
produced with a rate of about 7 % (Chiba, 1987) and kaons with a rate of about 6% of 
all annihilations (Sedlak and Simak, 1988). 

In fireball models the pion multiplicity N = N + + iV_ + Nq follows a Gaussian dis- 
tribution (Orfanidis and Rittenberg, 1973). The pion multiplicity distribution at rest in 
liquid hydrogen is shown in Fig. |l|. Following the model of Pais (1960) one expects on sta- 
tistical grounds the branching ratios to be distributed according to 1/(N + \N-\Nq\) for a 
given multiplicity N. The open squares show the predictions normalized to the measured 
branching ratios from channels with charged pions and Nq < 1 (Armenteros and French, 
1969). The full circles show the data from bubble chambers, together with Crystal Barrel 
results for iVo > 1. The fit to the data (curve) leads to a = 1 for a Gauss distribution 
assuming (N) = 5. The open circles show an estimate from bubble chamber experiments 
which appears to overestimate the contribution from N = 5 (Ghesquiere, 1974). 

2.1 S- and P-wave annihilation at rest 

Stopping antiprotons in hydrogen are captured to form antiprotonic hydrogen atoms (pro- 
tonium) . The probability of forming a pp atom is highest for states with principal quan- 
tum number n ~ 30 corresponding to the binding energy (13.6 eV) of the K-shell electron 
ejected during the capture process. Two competing de-excitation mechanisms occur: (i) 
the cascade to lower levels by X-ray or external Auger emission of electrons from neigh- 
bouring H2-molecules and (ii) Stark mixing between the various angular momentum states 
due to collisions with neighbouring H2 molecules. Details on the cascade process can be 
found in Batty (1989). In liquid hydrogen, Stark mixing dominates (Day, 1960) and the 
pp system annihilates with the angular momentum £ = from high S levels (S-wave an- 
nihilation) due to the absence of angular momentum barrier. The initial states are the 
spin singlet (s = 0) 1 Sq and the spin triplet (s = 1) 3 Si levels with parity P = (— 
and C-parity C = (— l)^ +s , hence with quantum numbers 

jPC^So) = 0~+ and J PC ( 3 5i) = 1—. (1) 

The cascade, important in low density hydrogen (e.g. in gas), populates mainly the 
n = 2 level (2P) from which the pp atom annihilates due its small size: The K a transition 
from 2P to IS has been observed at LEAR in gaseous hydrogen (Ahmad, 1985; Baker, 
1988). Compared to annihilation, it is suppressed with a probability of (98 ± 1) % 
at atmospheric pressure. Annihilation with relative angular momentum I = 1 (P-wave 
annihilation) can therefore be selected by detecting the L X-rays to the 2P levels, in 
coincidence with the annihilation products. This procedure permits the spectroscopy of 
intermediate meson resonances produced from the P-states 

J PC ( 1 P 1 ) = 1+-, J^( 3 P ) = ++ , 

J^(3 Pl ) = 1++ and J PC ( 3 P 2 ) = 2++. (2) 

Annihilation from P-states has led to the discovery of the /2(1565) meson (May, 1989). 

The much reduced Stark mixing in low density hydrogen also allows annihilation from 
higher P levels. At 1 bar, S- and P-waves each contribute about 50% to annihilation 
(Doser, 1988). Annihilation from D-waves is negligible due to the very small overlap of 
the p and p wave functions. 

The assumption of S-wave dominance in liquid hydrogen is often a crucial ingredient to 
the amplitude analyses when determining the spin and parity of an intermediate resonance 
in the annihilation process, since the quantum numbers of the initial state must be known. 
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The precise fraction of P-wave annihilation in liquid has been the subject of a longstanding 
controversy. The reaction pp — ► 7r°7r° can only proceed through the P-states ++ or 2 ++ 
(see section |4.1| ) while it + tt~ also proceeds from S-states (1 ). The annihilation rate 
B(ir°ir ) for this channel in liquid has been measured earlier by several groups but with 
inconsistent results (Devons, 1971; Adiels, 1987; Chiba, 1988). 

Crystal Barrel has determined the branching ratio for pp — ► 7r 7r in liquid by mea- 
suring the angles and energies of the four decay photons (Amsler, 1992a). The main 
difficulties in selecting this channel are annihilation into 3tt° which occurs with a much 
higher rate and, most importantly P-wave annihilation in flight. The former background 
source can be reduced with the good 7 detection efficiency and large solid angle of Crystal 
Barrel, while the latter can be eliminated thanks to the very narrow stop distribution from 
cooled low-energy antiprotons from LEAR (0.5 mm at 200 MeV/c). The small contami- 
nation from annihilation in flight can easily be subtracted from the stopping distribution 
by measuring the annihilation vertex. The latter was determined by performing a 5 con- 
straints (5C) fit to pp — > 7r°7T , assuming energy conservation, two invariant 27-masses 
consistent with 2tt and momentum conservation perpendicular to the beam axis. The 
branching ratio for 7r°7r° is 

B(tt°tt°) = (6.93 ± 0.43) x 10~ 4 , (3) 

in agreement with Devons (1971) but much larger than Adiels (1987) and Chiba (1988). 
From the annihilation rate B(ir + Tr~)2P into 7r + 7r~ from atomic 2P-states (Doser, 1988) 
one can, in principle, extract the fraction f p of P-wave annihilation in liquid 2 : 

f P = 2 „A- y = (28.8 ± 3.5)%. (4) 

B{lt+Tt ) 2 p 

This is a surprisingly large contribution. However, Eq. (Q) assumes that the population 
of the fine and hyperfine structure states is the same for the 2P as for the higher P 
levels which is in general not true. In liquid, strong Stark mixing constantly repopulates 
the levels. A P-state with large hadronic width, for instance 3 Po (Carbonell, 1989), will 
therefore contribute more to annihilation than expected from a pure statistical population. 
On the other hand, in low pressure gas or for states with low principal quantum numbers 
the levels are populated according to their statistical weights. The branching ratio for 
annihilation into a given final state is given in terms of the branching ratios Bf and Bf 
from the two S-, respectively the four P-states (Batty, 1996): 

B = {l-fp{p)]j^wfEf{p)Bf 

i=l 

+fp(p)j2w[E[(p)Bf, (5) 
i=i 

where p is the target density. The purely statistical weights are 

s 2Jj + 1 p 2Ji + 1 . . 

wf = , wf = . (6) 

4 12 

The enhancement factors Ei describe the departure from pure statistical population 
(Ei=l). For 7r°7T in liquid one obtains 

P(vrV) = /p(hq)[lp3 Po (liq)P3p (7rV) 



2 In Doser (1988) f p was found to be (8.6 ± 1.1) %, when using the most precise measurement for 7r°7r° 
available at that time (Adiels, 1987). 
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+^Es P2 (liq)Bs P2 ^\ )), (7) 

and for 7r + 7r _ from 2P states 

B{^-) 2P = 2[±B 3Po (tt\°) + ^Bs P2 (k\ )]. (8) 

It is obviously not possible to determine /p(liq) unless the enhancement factors are unity 
(and hence Eq. (^) follows). The enhancement factors have been calculated with an 
X-ray cascade calculation (Batty, 1996) using the observed yields of K and L X-rays in 
antiprotonic atoms and the predicted hadronic widths from optical potential models of the 
pp interaction (Carbonell, 1989). For example, Batty (1996) finds typically E3 Po (liq) ~ 
2.3 and i?3p 2 (liq) ~ 1.0. The branching ratios Bi and fp(p) were then fitted to the 
measured two-body branching ratios for pp — » 7r°7r°, 7r + 7r~, K + K~ , K$Ks and KsKl at 
various target densities, with and without L X-ray coincidence. The fraction of P-wave 
annihilation is shown in Fig. | as a function of density. In liquid hydrogen one obtains 

/ P (liq) = (13±4)%, (9) 

a more realistic value when compared to Eq. (Q). 



3 The Crystal Barrel Experiment 
3.1 Detector 

Figure |3| shows a sketch of the Crystal Barrel detector (Aker, 1992). The incoming 
antiprotons entered a 1.5 T solenoidal magnet along its axis and interacted in a liquid 
hydrogen target, 44 mm long and 17 mm in diameter. A segmented silicon counter in 
front of the target defined the incoming beam. The final state charge multiplicity was 
determined online with two cylindrical proportional wire chambers (PWC). The charged 
particle momentum was measured by a jet drift chamber (JDC) which also provided tt/K 
separation below 500 MeV/c by ionization sampling. 

Photons were detected in a barrel-shaped assembly of 1,380 CsI(Tl) crystals, 16.1 
radiation lengths long (30 cm), with photodiode readout. The crystals were oriented 
towards the interaction point and covered a solid angle of 0.97 x 47r. Each crystal, wrapped 
in teflon and aluminized mylar, was enclosed in a 100 ^m thick titanium container. The 
light (peaking at 550 nm) was collected at the rear end by a wavelength shifter and the 
re-emitted light was detected by a photodiode glued on the edge of the wavelength shifter. 
With the low electronic noise of typically 220 keV the energy resolution was 

; = do) 

and photons could be detected efficiently down to 4 MeV. The angular resolution was 
typically a = 20 mrad for both polar and azimuthal angles. The mass resolution was a 
= 10 MeV for tt° and 17 MeV for rj -> 2 7 . 

A rough calibration of the electromagnetic calorimeter was first obtained with travers- 
ing minimum ionizing pions which deposit 170 MeV in the crystals. The final calibration 
was achieved with 0-prong events using 27 invariant masses from 7r° decays. An energy 
dependent correction was applied to take shower leakage at the rear end of the crystals 
into account. The stability of the calibration was monitored with a light pulser system. 
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The JDC had 30 sectors, each with 23 sense wires at radial distances between 63 mm 
and 239 mm) read out on both ends by 100 MHz flash ADC's. The position resolution 
in the plane transverse to the beam axis (rcj) coordinates) was a = 125 /im using slow 
gas, a 90:10 % CC^/isobutane mixture. The coordinate z along the wire was determined 
by charge division with a resolution of a = 8 mm. This led to a momentum resolution 
for pions of a jp ~2 % at 200 MeV/c, rising to ~ 7% at 1 GeV/c for those tracks that 
traversed all JDC layers. 

The z coordinates were calibrated by fitting straight tracks from 4-prong events with- 
out magnetic field. The momentum calibration was performed with monoenergetic pions 
and kaons from the two-body final states 7r + 7r~ and K + K~. Pressure and temperature 
dependent drift time tables were generated and fitted to the measured momentum distri- 
bution. 

In 1994 the JDC was replaced by a new jet drift chamber with only 15 sectors for the 
6 innermost layers. In 1995 the PWC's were also replaced by a microstrip vertex detector 
(SVX) consisting of 15 single-sided silicon detectors arranged in a windmill configuration 
at a radial distance of 13 mm around the target (Fig. ||). Each detector had 128 strips with 
a pitch of 50 ^m running parallel to the beam axis. The increase of charge multiplicity 
between the SVX and the inner layers of the JDC permitted to trigger on K$ — ► ir + ir~ . 
The SVX also provided an improved vertex resolution in r<j) and a better momentum 
resolution. 

For annihilation at rest in liquid hydrogen the p incident momentum was 200 MeV/c 
with typically 10 4 incident p/s to minimize pile-up in the crystals. For annihilation in 
gaseous hydrogen the liquid target was replaced by a hydrogen flask at 13 bar. The 
incident momentum was 105 MeV/c. Since the annihilation rate was higher than the 
maximum possible data acquisition speed, a multilevel trigger could be used. The two 
PWC's and the inner layers (2 - 5) of the JDC determined the charged multiplicity of the 
final state. Events with long tracks could be selected for optimum momentum resolution 
by counting the charged multiplicity in the outer layers (20 and 21) of the JDC. A hard- 
wired processor determined the cluster multiplicity in the barrel. A processor then fetched 
the digitized energy deposits in the barrel, computed all two-photon invariant masses thus 
providing a trigger on the 7r° or rj multiplicity (Urner, 1995). 

3.2 Photon reconstruction 

We now briefly describe the photon reconstruction which is particularly relevant to the 
results reviewed in this article. Photon induced electromagnetic showers spread out over 
several crystals. The size of a cluster depends on the photon energy and varies from 1 
to about 20 crystals. The reconstruction of photons is done by searching for clusters of 
neighbouring crystals with energy deposits of at least 1 MeV. The threshold for cluster 
identification (typically between 4 and 20 MeV) depends on the annihilation channel being 
studied. Local maxima with a predefined threshold (typically between 10 and 20 MeV) 
are then searched for within clusters. When only one local maximum is found, the photon 
energy is defined as the cluster energy and the direction is given by the center of gravity of 
the crystals, weighted by their energies. When n local maxima are found within a cluster, 
the latter is assumed to contain showers from n photons. In this case the cluster energy 
Eq is shared between the n subclusters of nine crystals with energies E% around the local 
maxima. Hence the photon energies are given by 

E Jti = ^-etE c . (11) 
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Additional clusters mocking photons are due to shower fluctuations which may develop 
small but well separated satellites in the vicinity of the main shower. These "split-offs" 
can be removed by requiring a minimum separation between the showers. However, this 
cut may reduce the detection efficiency for high energy 7r°'s since photons from ir° decay 
cluster around the minimum opening angle. The opening angle between two photons with 
energies E\ < E2 from tt° decay is given by 

(l + R) 2 E 1 , s 

cos0 = l-i^f withi?=-^, (12) 

where 7 = E n o / m w o . Hence for all pairs of neighbouring clusters one calculates R and 
removes the low energy clusters whenever cos</> is larger than given by Eq. (|T2|), assuming 
the maximum possible value of 7 in the annihilation channel under consideration (Pietra, 
1996). 

Clusters generated by ionizing particles can be removed by matching the impact points 
extrapolated from the reconstructed tracks in the JDC. However, split-offs from charged 
particles are more cumbersome to eliminate. They are initiated, for example, from neu- 
trons which travel long distances before being absorbed. These split-offs can be suppressed 
by requiring momentum and energy conservation in the annihilation process (kinematic 
fits). 



3.3 Available data 

The bulk of the Crystal Barrel data consists in pp annihilation at rest and in flight in 
liquid hydrogen. As discussed above, annihilation from initial P-states is enhanced when 
using a gaseous target. Annihilation in deuterium at rest allows the formation of NN 
bound states below 2mjv> the spectator neutron (or proton) removing the excess energy 
(for a review on baryonium states, see Amsler (1987)). With a spectator proton one gains 
access to pn annihilation, a pure isospin 1 = 1 initial state. 

The data collected by Crystal Barrel are shown in Table Data were taken in liquid 
hydrogen, gaseous hydrogen (13 bar) and in liquid deuterium with a minimum bias trigger 
(requiring only an incident antiproton) or with the multiplicity trigger requiring 0-prong 
or n-prong with long tracks in the JDC. In addition, data were collected with specialized 
triggers enhancing specific final states. As a comparison, the largest earlier sample of 
annihilations at rest in liquid was obtained by the CERN-College de France collaboration 
with about 100,000 pionic events and 80,000 events containing at least one K$ — ► ir + iT~ 
(Armenteros and French, 1969). The Asterix collaboration collected some 10 7 pionic 
events in gaseous hydrogen at 1 bar (for a review and references see Amsler and Myhrer 
(1991)). The total number of annihilations at rest in liquid hydrogen collected by Crystal 
Barrel is 10 8 . The triggered 0-prong sample alone corresponds to 6.3xl0 8 annihilations. 



4 Annihilation into Two Mesons 

Consider a pair MM of charge conjugated mesons in the eigenstate of isospin /. The P-, 
C- and G-parities are: 

P(MM) = (-1) L , (13) 
C(MM) = (-1) L+S , (14) 
G{MM) = (-1) L+S+I , (15) 
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where L is the relative angular momentum and S the total spin. For the pp system with 
angular momentum £ and spin s one has 

p(pp) = (-r/+\ (i6) 

cm = (-i) e+s , (17) 

G(pp) = (-lY +s+I . (18) 

For annihilation into two mesons the two sets of equations relate the quantum numbers of 
the initial state to those of the final state since P, C, G and I are conserved. In addition, 
L, S, I and s must be chosen so that the total angular momentum J is conserved: 

\L-S\<J<L + S, \l-s\<J<£ + s. (19) 

Since P, C and G are multiplicative quantum numbers these relations are especially 
restrictive for mesons that are eigenstates of C and G, e.g. for neutral non-strange mesons. 
For example, for two identical neutral non-strange pseudoscalars (e.g. 7r 7r°) with S = 0, 
C = +1, Eq. ( |l4|) implies that L is even and then Eq. (li~6| ) requires £ to be odd 
(annihilation from P-states only). Equation (pi]) further requires with G = +1 that 1 = 
and hence with Eq. ( |l8| ) annihilation from the (/ = 0) ++ or 2 ++ atomic states. 

For a pair of non-identical neutral pseudoscalars (e.g. 7t°t/) L may be odd and hence 
the possible quantum numbers are ++ , 1 2 ++ , 3 etc. However, 1 h and 3 h do 
not couple to pp since Eqs. ([Hf ) and (|l7]) require i even and s = and hence J even. In 
fact these "exotic" quantum numbers do not couple to any fermion-antifermion pair and 
are, in particular, excluded for qq mesons. 



4.1 Annihilation into two neutral mesons 

Crystal Barrel has measured the branching ratios for pp annihilation into two neutral light 
mesons from about 10 7 annihilations into 0-prong (Amsler, 1993b). These data have been 
collected by vetoing charged particles with the PWC's and the internal layers of the JDC. 
The lowest 7-multiplicity was four (e.g. 77) and the highest nine (e.g. -quj, with 

77 — > 37T° and u> — ► 7r°7). To control systematic errors in the detection efficiency, some of 
these branching ratios have been determined from different final state multiplicities. For 
example, 77 decays to 27 and 3ir° and hence r\r\ is accessible from 47 and 87 events. 

Figure [5] shows a scatterplot of 27- invariant masses for events with 47. Events have 
been selected by requiring four clusters in the barrel and applying momentum and energy 
conservation (4C fit). Signals from 7r°7r , vr°r/, 7777 and even 7r ?/ are clearly visible. The 
dark diagonal band at the edge is due to wrong combinations. The detection and recon- 
struction efficiency was typically 40% for 47 events, obtained by Monte Carlo simulation 
with GEANT. As discussed above, two neutral pseudoscalars couple only to atomic P- 
states and are therefore suppressed in liquid hydrogen. On the other hand, the channels 
and r/co couple to 3 Si and hence have a larger branching ratio. In spite of the good 
detection efficiency of the detector one therefore observes the background signals from 
and rjuj, where oj decays to with a missing (undetected) photon. 

Figure |6| shows the 7r°7 momentum distribution for pp — > 47r°7 events (8C fit requiring 
47T°). The peak at 657 MeV/c is due to the channel pp — ► rjuj. For these 97 events the 
detection efficiency was 10%. 

The branching ratios are given in Table |2|. They are always corrected for the unob- 
served (but known) decay modes of the final state mesons (Barnett, 1996). For Crystal 
Barrel data the absolute normalization was provided by comparison with 7r°7r° which has 
been measured with minimum bias data (Eq. (§)). 
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Signals for uj-k and ur) have also been observed for to decaying to ir + ir~ir°, leading 
to 7r + 7r~47 (Schmid, 1991). Figure ^ shows the tt + tt^tt° invariant mass spectrum for 
7r + 7r _ 7r°ry events. The branching ratio for ujtt and ior\ are in excellent agreement with the 
ones from 0-prong (Table [2]). 

The angular distribution in the uo rest frame contains information on the initial atomic 
state. The distribution of the angle between the normal to the plane spanned by the three 
pions and the direction of the recoiling r\ is plotted in the inset of Fig. 0. Using the method 
described in section 7.2 one predicts the distribution sin 2 # for annihilation from 3 SI while 
the distribution should be isotropic for annihilation from P\. The fit (curve) allows (12 
±4)% P-wave. Figure || shows the angular distribution of the 7 in the u rest frame for 
lut}(lu — > 7r°7). The predicted distribution is (1 + cos 2 6) from 3 Si and is again isotropic 
for 1 P i . The fit (curve) allows (9 ± 3)% P-wave. However, these results assume that the 
relative angular momentum between r\ and uj is L = from Pi, thus neglecting L = 2. 
Without this assumption, the fraction of P-wave cannot be determined from the angular 
distributions due to the unknown interference between the L = and L = 2 amplitudes. 

Some of the branching ratios for two-neutral mesons have been measured earlier 
(Adiels, 1989; Chiba, 1988) by detecting and reconstructing 7r°'s or rj's with small solid 
angle detectors and observing peaks in the 7r° or 77 inclusive momentum spectra. Since 
the branching ratios are small, these early data are often statistically weak or subject 
to uncertainties in the baseline subtraction from the inclusive spectra. In fact most of 
the Crystal Barrel results disagree with these measurements which should not be used 
anymore. Table therefore updates Table 1 in Amsler and Myhrer (1991). 



4.2 The annihilation mechanism 

There is currently no model which completely and satisfactorily describes the measured 
two-body branching ratios listed in Table ^ (for a review of annihilation models and refer- 
ences, see Amsler and Myhrer (1991). Since the proton and the antiproton wavefunctions 
overlap one expects quarks to play an important role in the annihilation dynamics. For 
instance, pp annihilation into two mesons can be described by the annihilation of two qq 
pairs and the creation of a new pair (annihilation graph A) or by the annihilation of one 
qq pair and the rearrangement of the other two pairs (rearrangement graph R) , see Fig. ||. 
At low energies there is, however, no consensus as to which operator should be used to 
describe the emission and absorption of gluons. In a first approach one assumes that only 
the flavor flow between initial and final states is important (Genz, 1983; Hartmann, 1988). 
The Quark Line Rule (QLR) states that annihilation into uu and dd is excluded if the A 
graph dominates, while annihilation into two dd mesons is forbidden if R dominates (see 
Fig. |). The OZI-rule (Okubo, 1963) is a special case of the QLR: Annihilation into one 
or more ss mesons is forbidden. We shall confront these simple rules below with Crystal 
Barrel data. 

Another approach, which we shall use, is the nearest threshold dominance model which 
describes reasonably well the observed final state multiplicity as a function of p momentum 
(Vandermeulen, 1988). The branching ratio for annihilation into two mesons with masses 
m a und 77ib is given by 

W = P C C ab exp(-A v /s- (m a + m b ) 2 ), (20) 

where p is the meson momentum in the pp center of mass system with total energy ^/s, 
Co a normalization constant and C a b a multiplicity factor depending on spin and isospin. 
The constant A= 1.2 GeV -1 has been fitted to the cross section for pp annihilation into 
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as a function of p momentum. For annihilation into kaons the fit to kaonic channels 
requires the additional normalization factor C\/Cq = 0.15. Thus annihilation into the 
heaviest possible meson pair is enhanced with respect to phase space p by the exponential 
form factor in Eq. (|20|). This is natural in the framework of baryon exchange models 
which prefer small momentum transfers at the baryon-meson vertices. 

In more refined models the branching ratios for annihilation at rest depend on the 
atomic wave function distorted by strong interaction at short distances (Carbonell, 1989). 
Predictions for the branching ratios therefore depend on models for the meson exchange 
potential which are uncertain below 1 fm. Also, the quark description has to be comple- 
mented by baryon and meson exchanges to take the finite size of the emitted mesons into 
account. 

In the absence of strong interaction the pp atomic state is an equal superposition of 
isospin 1 = and 1 states. Naively one would therefore expect half the protonium states 
to annihilate into a final state of given isospin. However, pp to nn transitions occuring 
at short distances may modify the population of I = and 1 = 1 states and therefore 
enhance or reduce the annihilation rate to a final state of given isospin (Klempt, 1990; 
Jaenicke, 1991). Nonetheless, one expects that predictions for ratios of branching ratios 
for channels with the same isospin and proceeding from the same atomic states are less 
sensitive to model dependence. We shall therefore compare predictions from the QLR 
with ratios of branching ratios from Table [|. 

The flavor content of the rj and rj' mesons is given by 



= —J={\ uu ) + \dd)) sin(0j - 9 p ) - \ss) cos(#j - 9 P ), 
a/2 

W) = ^(\uu) + \dd))cos(9 i -9 p ) + \ss)sm{9 i -9 p ), (21) 



where 9i = 35.3° is the ideal mixing angle. The flavor wave functions of the 7r° and p° are 

|A|p°) = -^(|^)-M), (22) 
and those of u and 4>, assuming ideal mixing in the vector nonet, 

\u) = ^=(\uu) + \dd)), \<f>) = -\ss). (23) 

The branching ratio for annihilation into two neutral mesons is then given by B = B ■ W 
with 

B = \{pp\T\M l M 2 )\ 2 = iXJrCMi, te^] 2 )te|M 1 )(^|M 2 )| 2 , (24) 

and | qi = u or d (Genz, 1985). In the absence of ss pairs in the nucleon, the QLR 
forbids the production of ss mesons and therefore the ss components in Mi and M 2 can 
be ignored in Eq. (|24|). The predicted ratios of branching ratios are given by the first 
four rows in Table ^ for various channels from the same atomic states. They depend only 
on the pseudoscalar mixing angle 9 P . To extract 9 p the measured branching ratios from 
Table | must be first divided by W (Eq. @), ignoring Cq and C a b which cancel in the 
ratio. 

The pseudoscalar mixing angle has been measured in various meson decays (e.g. r\ and 
rj' radiative decays, J/tp radiative decays to rj and rj') and is known to be close to -20° 



3 The theoretical prediction B has to be multiplied by two for a pair of non-identical mesons. 
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(Gilman and Kauffmann, 1987). The agreement with our simple model of annihilation is 
amazing (third column of Table ^) . We emphasize that the predictions in the upper four 
rows of Table || are valid independently of the relative contributions from the A and R 
graphs. 

Conversely, one can assume the validity of the model and extract from the first four 
rows in Table ^ the average 

p = (-19.4 ±0.9)°. (25) 

Leaving the constant A in Eq. (^) as a free fit parameter one obtains 9 P = (— 17.3 ± 1.8)° 
from early Crystal Barrel data (Amsler, 1992b). Assuming now dominance of the planar 
graph A, the amplitudes T([uH]i, [fid^) and T([dd]i, [uu]2) vanish and one obtains the 
predictions in the lower part of Table ||. The measurements lead in general to incorrect 
values for 6 P , presumably due to the contribution of the R graph. Also, for p°p° and lulu 
one expects from A dominance 

B(p°p°) = B{lulu), (26) 

in violent disagreement with data (Table ^). On the other hand, if R dominates the 
amplitude T([dd]\ : [dd]2) vanishes and one predicts from Eq. ( p4[ ) with a 2 = sin 2 (#j — 9 P ) 
the inequality (Genz, 1990) 



|a 2 y2B(7r°7r°) - ^2B(rjrj)\ 2 < Aa B(tt rj) < \a 2 : y '2B(tt°tt ) + \J2B(r]rj)\ 2 , (27) 

which is fulfilled by data. 

There is, however, a caveat: the predictions ( pl|) have been compared to the measured 
branching ratios corrected by W. As pointed out earlier, Eq. ([H]) provides a good fit to 
the mutiplicity distribution in low energy pp annihilation as a function of p momentum. 
Other correcting factors can, however, be found in the literature. In section |7| we shall 
use the phase space factor 

W = pFl(p) (28) 

where Fi(p) is the Blatt-Weisskopf damping factor which suppresses high angular mo- 
menta L for small p. This factor is determined by the range of the interaction, usually 
chosen as 1 fm (pr =197 MeV/c). Convenient expressions for Fl(p) are given in Table ||. 
For p much larger than pr, Fi(p) — 1 and for p much smaller than pr 

F L (p)~p L . (29) 

This last prescription provides a reasonable agreement when comparing the measured 
decay branching ratios of mesons, especially tensors, with predictions from SU(3), as we 
shall discuss in section [10.1 , These alternative phase space factors may also be used to 



determine the pseudoscalar mixing angle. However, they do not lead to consistent values 
for 6 P (Amsler, 1992b). Agreement is achieved with prescription (p0[), which we shall also 
employ in the next section. 

In conclusion, the naive quark model assuming only the QLR and two-body threshold 
dominance reproduces the correct pseudoscalar mixing angle from the measured two- 
meson final states. This is a clear indication for quark dynamics in the annihilation 
process. The relative contribution from R and A cannot be extracted but the non-planar 
graph R must contribute substantially to the annihilation process. 
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5 Electromagnetic Processes 



5.1 Radiative annihilation 

For annihilations leading to direct photons the isospin 1 = and 1 amplitudes from the 
same pp atomic state interfere since isospin is not conserved in electromagnetic processes. 
Radiative annihilation pp — > jX, where X stands for any neutral meson, involves the 
annihilation of a qq pair into a photon. The branching ratios can be calculated from the 
Vector Dominance Model (VDM) which relates 7 emission to the emission of p, lo and (ft 
mesons (Delcourt, 1984). Assuming ideal mixing in the vector nonet one may actually 
neglect (ft production which is forbidden by the QLR. The amplitude for pp — > ^X is then 
given by the coherent sum of the two 1 = and 1 = 1 amplitudes with unknown relative 
phase j3 (fig. |l^). According to VDM, the -jp coupling g p7 is a factor of three stronger 
than the ju> coupling. The branching ratio is then^ 



with 



B(rfX) = A 2 [B(pX) + ±B{uX) + B{pX)B{uX) cos/3], (30) 



A = 2E£1 = 0.055. (31) 



Equation ( |30| ) provides lower and upper limits (cos/3 = ± 1) for branching ratios. 

Radiative annihilation has not been observed so far with the exception of 7r°7 (Adiels, 
1987). Crystal Barrel has measured the rates for 7r°7, 777, W7 and has obtained upper 
limits for 7/7 and 77 (Amsler, 1993c). Annihilation into 07 (Amsler, 1995c) is treated in 



section 6.1 



The starting data sample consisted of 4.5 x 10 0-prong events. Figure 11 shows a 
typical 7r°7 event leading to three detected photons. The main background to 7r°7 stems 
from annihilations into vr vr for which one of the photons from 7r° decay has not been 
detected, mainly because its energy lies below detection threshold (10 MeV). The 7r° 
momentum for 7T°7 is slightly higher (5 MeV/c) than for 7r°7r°. The small downward shift 
of the 7T° momentum peak due to the contamination could be observed thanks to 
the good energy resolution of the detector and could be used to estimate the feedthrough 
from vr vr : (29 ± 8) %, in agreement with Monte Carlo simulations. The result for the 
7r°7 branching ratio is given in Table |. It disa grees with the one obtained earlier from 
the 7T° inclusive momentum spectrum: (1.74 ±0.22) x 10 -5 (Adiels, 1987). 

Results for 777,0^7 and 7/7 are also given in Table |. The uj was detected in its 7r°7 
and the rf searched for in its 27 decay mode. The main contaminants were r?7r°, lott and 
t]'tt°, respectively, with one photon escaping detection. For pp — » 77, 98 ± 10 events were 
observed of which 70 ± 8 were expected feedthrough from 7r°7 and 7r°7r°. This corresponds 
to a branching ratio of (3.3 ±1.5) x 10 -7 which the collaboration prefers to quote as an 
upper limit (Table ||). 

The branching ratios, divided by W (Eq. (^(J)), are compared in Table |5| with the 
range allowed by Eq. (|30|). Apart from <ftj to which we shall return later, the results 
agree with predictions from VDM. For and 777 (from 3 Si) the isospin amplitudes 
interfere destructively (cos/3 ~ -0.3) while for loj (from x 5o) they interfere constructively 
(cos/3 ~ 0.13). We emphasize that these conclusions depend on the prescription for the 
phase space correction. With a phase space factor of the form p 3 one finds strongly 
destructive amplitudes (Amsler, 1993c), see also Locher (1994) and Markushin (1997). 



4 For two identical particles, e.g. p°p° or uilu, the measured branching ratios divided by phase space, B, have 
to be multiplied by two. 
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No prediction can be made from VDM for 77 due to the contribution of three ampli- 
tudes with unknown relative phases: p°uj from 1=1, p°p° and loco from 1 = 0. Also, the 
branching ratio for p°p° is poorly known (Table ||). 

5.2 Search for light gauge bosons in pseudoscalar meson 
decays 

Extensions of the standard model allow additional gauge bosons, some of which could 
be light enough to be produced in the decay of pseudoscalar mesons (Dobroliubov and 
Ignatiev, 1988). Radiative decays 7T°, 77, 7/ — ► 7X are particularly suitable since they are 
only sensitive to gauge bosons X with quantum numbers J p = 1~ . Branching ratios are 
predicted to lie in the range 10 -7 to 10 -3 (Dobroliubov, 1990). Experimental upper limits 
for 7T° — > jX are of the order 5 x 10~ 4 for long lived gauge bosons with lifetime r > 10 
s (Atiya, 1992). Short lived gauge bosons decaying subsequently to e + e~ are not observed 
with an upper limit of 4 x 10~ 6 (Meijer Drees, 1992). 

Crystal Barrel has searched for radiative decays where X is a long lived weakly inter- 
acting gauge boson escaping from the detector without interaction, or decaying to vv. The 
search was performed using the reactions pp — ► 3ir°, 7r°7r°r/ and 7r 7r°r/' at rest (Amsler, 
1994a, 1996a) which occur with a sufficiently high probability (see Table ^ below) and 
are kinematically well constrained. Events with five photons were selected from a sample 
of 15 million annihilations into neutral final states (18 million for rj' decays). Since the 
branching ratio for 0-prong annihilation is about 4%, the data sample corresponds to some 
400 million pp annihilations in liquid hydrogen. 

Events consistent with decays and a single (unpaired) 7 were then selected by 

requiring energy and momentum conservation for pp annihilation into 7r°7r 7r°, 7r°ir°r] or 
7r°7ry with a missing 7r°, 77 or rj . Thus 3C kinematic fits were applied, ignoring the 
remaining fifth photon. The measured energy of the latter was then transformed into the 
rest frame of the missing pseudoscalar. In this frame a missing state X with mass mx 
would produce a peak in the 7-energy distribution at 

(32) 

with width determined by the experimental resolution, where m is the mass of the missing 
pseudoscalar. Thus, if X is simply a missing (undetected) 7 from 7r°,7/ or rj' decay, one 
finds with mx = that E* = m/2, as expected. 

The main source of background is annihilation into three pseudoscalars for which one of 
the photons escaped detection. This occurs for (i) photons with energies below detection 
threshold {E < 20 MeV) or (ii) for photons emitted into insensitive areas of the detector. 
The latter background can be reduced by rejecting events for which the missing 7 could 
have been emitted e.g. in the holes along the beam pipe. The high efficiency and large 
angular coverage of the Crystal Barrel are therefore crucial in this analysis. For the 3n° 
channel an important background also arose from pp — > tt°lo with uj — > 7r°7, leading to 
57. This background could be reduced by rejecting events which satisfy the 7r°u;(— > 7r°7) 
kinematics. For the background channel K${-^ 2it°)Kl with an interacting Kl 

faking a missing 77 could be eliminated with appropriate kinematic cuts. 

The -E*-energy distribution for 7r° decay is shown in Fig. 12. The broad peak 
around 70 MeV is due to 7r° decays into 27 where one 7 has escaped detection, or from 
residual 7r°u; events. The fit to the distribution (full line) agrees with the simulated rate 
of background from 3tt° and The dotted line shows the expected signal for a state 
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with mass mx = 120 MeV, produced in tt° decay with a branching ratio of 5 x 10 -4 . The 
corresponding distributions for rj and rf decays can be found in Amsler (1996a). 

Upper limits for radiative decays are given in Fig. 13 as a function of mx- The upper 
limit for 7r° decay is an order of magnitude lower than from previous experiments (Atiya, 
1992). For r\ and rf decays no limits were available previously. Light gauge bosons are 
therefore not observed in radiative pseudoscalar decays at a level of 10 -4 to 10 -5 . 



5.3 r; — > 37r 

The 3ir decay of the r/ plays an important role in testing low energy QCD predictions. 
This isospin breaking decay is mainly due to the mass difference between u and d quarks. 
Crystal Barrel has measured the relative branching ratios for rj — ► 7r + 7r~7r°, rj —* 3tt° and 
rj — > 27 from samples of annihilations into 27r + 27r _ 7r°, 7r + 7r~37r° and 7r + 7r~27, respectively 
(Amsler, 1995d). The ratios of partial widths are 

n - tP^% = 1-44 ±0.13, 

1 (77 — > 7T+7T 7T U ) 

r 2 = r( W t _ , =1-78 ±0.16. (33) 

r(77 — > 7T+7T 7T U ) 

The result for n is in good agreement with chiral perturbation theory: 1.43 db 0.03 
(Gasser and Leutwyler, 1985) and 1.40 ± 0.03 when taking unitarity corrections into 
account (Kambor, 1996). With the known 27 partial width (Barnett, 1996) one can 
calculate from T2 the partial width r(7r + 7r~7r°) = 258 db 32 eV, in accord with chiral 
perturbation theory (r = 230 eV), taking into account corrections to the u — d mass 
difference (Donoghue, 1992). In good approximation, the rj — ► 7r + 7r~7r° Dalitz plot may 
be described by the matrix element squared 

\M(x,y)\ 2 oc 1 + ay + by 2 (34) 

myq) — m(7r u ) — 2m(7r ± ) 
where Tq is the kinetic energy of the neutral pion. The parameters a and b were determined 
in Amsler (1995d), but more accurate values are now available from the annihilation 
channel pp — ► 7r°7r°r/. Abele (1997b) finds 

a = -1.19 ±0.07, b = 0.19 ±0.11, (36) 

in reasonable agreement with chiral perturbation calculations which predict a = -1.3 and 
b = 0.38 (Gasser and Leutwyler, 1985). 

The matrix element for rj decay to 3ir° is directly connected to the matrix element for 
the charged mode because 37r° is an / = 1 state. The matrix element squared for r\ decay 
to 3-7T is given by 

\M(z)\ 2 = l + 2az, (37) 
where z is the distance from the center of the rj — > 3tt° Dalitz plot, 



1=1 



. \ 1 1 — 1 1 ) 1 1 1 1 

(38) 



2 ^ f 3JS?i — m(r?) 



7/4(77) — 3m(7T ) 



and where E{ are the total energies of the pions. Chiral perturbation theory up to next- 
to-leading order predict a to be zero (Gasser and Leutwyler, 1985) leading to a homoge- 
neously populated Dalitz plot. Taking unitarity corrections into account, Kambor (1996) 
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predicts a ~ —0.01. Experiments have so far reported values for a compatible with zero, 
e.g. Aide (1984) finds -0.022 ± 0.023. 

Crystal Barrel has analyzed the 37r° Dalitz plot with 98,000 rj decays from the anni- 
hilation channel -K ir°r], leading to 10 detected photons (Abele, 1997c). The background 
and acceptance corrected matrix element is shown in Fig. [l4] as a function of z. The slope 
is clearly negative: 

a = -0.052 ± 0.020. (39) 



5.4 rf — > 7T + 7T 7 

The 7r + 7r~7 decay mode of the rf is generally believed to proceed through the p(770)7 
intermediate state (Barnett, 1996). However, the p mass extracted from a fit to the 
7r + 7r~ mass spectrum appears to lie some 20 MeV higher than for p production in e + e _ 
annihilations. This effect is due to the contribution of the direct decay into 7r + 7r~7 
(Bityukov, 1991) through the so- called box anomaly (Benayoun, 1993). Crystal Barrel 
has studied the rf — > 7r + 7r _ 7 channel where the rf is produced from the annihilation 
channels 7r°7r r/', ir + Tr~r]' and ujrf (Abele, 1997i). Evidence for the direct decay was 
confirmed at the 4<r level by fitting the 7r + 7r~ mass spectrum from a sample of 7,392 
rf decays. Including contributions from the box anomaly, the p mass turns out to be 
consistent with the standard value from e + e~ annihilation. Using the known two-photon 
decay widths of 77 and rf and the rj — ► 7r + 7r~7 decay spectrum from Layter (1973) the 
collaboration derived the pseudoscalar nonet parameters f n / fi=0.91 ± 0.02, /§/ '/ 7r =0.90 
± 0.05 and the pseudoscalar mixing angle 9 p = (—16.44 ± 1.20)°. 



5.5 Radiative uj decays 

The rates for radiative meson decays can be calculated from the naive quark model using 
SU(3) and the OZI rule (O'Donnell, 1981). Assuming ideal mixing in the vector nonet 
one finds, neglecting the small difference between u and d quark masses: 

B(u -» rn) _ 1 Pv 2 



cos^(54.7° + 6 P ) = 0.010, (40) 



where p n = 379 MeV/c and p v = 199 MeV/c are the decay momenta in the u rest frame 
and 9p is the pseudoscalar mixing angle (Amsler, 1992b). However, the production and 
decay of the to and p mesons are coupled by the isospin breaking uj to p transition since 
these mesons overlap (for references on p — uj mixing, see O'Connell (1995)). Since the 
width of p is much larger than the width of uj, the effect of p — uj mixing is essential in 
processes where p production is larger than uj production, for example in e + e~ annihilation 
where u and p are produced with a relative rate of 1/9. The determination of the branching 
ratio for uj — > 777 varies by a factor of five depending on whether the interference between 
w V7 an d P° — ► f77 is constructive or destructive (Dolinsky, 1989). A similar effect is 
observed in photoproduction (Andrews, 1977). The GAMS collaboration has determined 
the uj — > 777 decay branching ratio, (8.3 ± 2.1) x 10~ 4 , using the reaction ir~p — ► ujn at 
large momentum transfers, thus suppressing p production (Aide, 1994). 

In pp annihilation at rest the branching ratio for ujtt is much smaller than the branch- 
ing ratio for p°7r° while the converse is true for rjuj and r/p° (Table |||). We therefore 
expect that a determination of the uj — > 777 branching ratio from 7r° 'uj, neglecting p — uj 
mixing, will lead to a larger value than from 77a;. However, a simultaneous analysis of 
both branching ratios, including p — uj mixing, should lead to consistent results and allow 
a determination of the relative phase between the two amplitudes. 
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In Abele (1997d) the channels it uj and rjui (to —* 777) were reconstructed from 15.5 
million 0-prong events with five detected 7's. The events were submitted to a 6C kinematic 
fit assuming total energy and total momentum conservation, at least one 7r° — > 27 (or one 
77 — > 27) and to —* 37. The uj —* 37 Dalitz plots for it uj and 770; are shown in Fig. [D| 
using the variables 

T 2 -T, T 3 1 

x = ^7!^' y= Q-3' (41) 

where T±,T2 and T3 are the kinetic energies of the 7's in the uj rest frame and Q = 
T\ +T2 + I3. The prominent bands along the boundaries are due to w — > tt°j and the 
weaker bands around the center to uj — > r/j. 

The main background contributions arose from 67 events (pp — ► 37r°, 2ir°7], 2r]ir°) with 
a missing photon. This background (10% in the tt°uj and 5% in the 770; Dalitz plots) 
was simulated using the 3-pseudoscalar distributions discussed in section [8| and could be 
reduced with appropriate cuts (Pietra, 1996). After removal of the 7r° bands 147 ± 25 
w ~~ > VI events were found in the 77 bands of Fig. |l5|(a) and 123 ± 19 events in the 77 
bands of Fig. gj(b). 

The branching ratio for u — > 777 was derived by normalizing on the known branching 
ratio for u — > vr°7, (8.5 ± 0.5) % (Barnett, 1996). Correcting for the reconstruction 
efficiency Abele (1997d) finds a branching ratio of (13.1 ± 2.4) x 10~ 4 from and 
(6.5 ± 1.1) x 10~ 4 from ijlo, hence a much larger signal from tt°lo. 

Consider now the isospin breaking electromagnetic p — iv transition. The amplitude S 
for the reaction pp — » X(p — uj) — > A7/7 is, up to an arbitrary phase factor (Goldhaber, 
1969): 

S = 




+ ei(a+0) r^p ^_|^Lr_rj, (42) 

where ^4 is the production and T the decay amplitude of the two mesons and P p = 
m—m p +iT p /2, = The parameter 5 was determined from ui, p — > 7r + 7r _ : 

<5 = (2.48 ± 0.17) MeV (Weidenauer, 1993). The relative phase between the production 
amplitudes A p and A w is a while the relative phase between the decay amplitudes T p and 
is (p. In the absence of p — uj interference (5 = 0) Eq. ( p2[ ) reduces to a sum of two 
Breit-Wigner functions with relative phase a + <p. The magnitudes of the amplitudes A 
and T are proportional to the production branching ratios and the partial decay widths, 
respectively. 

The production phase a can be determined from p, uj — > ir + ir~ since the isospin vio- 
lating decay amplitude T(uj — > 7r + 7r~) may be neglected, leaving only the first term in Eq. 
(|42|). A value for a consistent with zero, (— 5.4±4.3)°, was measured by Crystal Barrel, us- 
ing the channel pp — ► t]tt + tt~ where p — uj interference is observed directly (Abele, 1997a). 
This phase is indeed predicted to be zero in e + e~ annihilation, in photoproduction and 
also in pp annihilation (Achasov and Shestakov, 1978). 

With the branching ratios for ujtt , ujt], ir°p° and rjp° given in Table |2| the intensity 
\S\ 2 was fitted to the number of observed u; — > r/7 events in ujtt and ujt], using Monte 
Carlo simulation. Both channels lead to consistent results for 

B(lo 777) = (6.6 ± 1.7) x 10~ 4 , (43) 

in agreement with the branching ratio from rjuj, obtained by neglecting p — uj interference. 
The phase 4> = (— 20^59)° l ea ds to constructive interference. The result Eq. (f43l) is in 
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excellent agreement with Aide (1994) and with the constructive interference solution in 
e + e~, (7.3 ±2.9) x (Dolinsky, 1989). This then solves the longstanding ambiguity in 
e + e~ annihilation between the constructive ((ft = 0) and destructive ((ft = ir) interference 
solutions. The branching ratio for p° — > 777, (12.2 ± 10.6) x 10 -4 , is not competitive but 
agrees with results from e + e~, (3.8 ± 0.7) x 10 -4 for constructive interference (Dolinsky, 
1989; Andrews, 1977). Using Eq. @ one then finds 

B ^^Vl) _ (7 . 8±2 . 1)xl0 -3 5 (44) 



B(u -» 7T°7) 

in agreement with SU(3). 

The u> — > 37 Dalitz plot is also useful to search for the direct process oj — > 37 which 
is similar to the decay of ( 3 S*i) orthopositronium into 37 and has not been observed so 
far. By analogy, the population in the u — > 37 Dalitz plot is expected to be almost 
homogeneous except for a slight increase close to its boundaries (Ore and Powell, 1949). 
Using the central region in Fig. |D](a) which contains only one event (6 entries) one obtains 
the upper limit B(uj — > tt°j) = 1.9 x 10~ 4 at 95% confidence level. This is somewhat more 
precise that the previous upper limit, 2 x 10 -4 at 90% confidence level (Prokoshkin and 
Samoilenko, 1995). 



6 Production of Mesons 

It has been known for some time that (ft production is enhanced beyond expectation from 



the OZI rule in various hadronic reactions (Cooper, 1978). Let us return to Eq. (21) 
and replace 77 by (ft and 7/ by u. The mixing angle becomes the mixing angle 9 V in the 
vector nonet. According to the OZI rule, (ft and uj can only be produced through their 
uu + dd components. Hence (ft production should vanish for an ideally mixed vector nonet 
(9 V = 0i) in which (ft is purely ss. Since (ft also decays to 3ir this is not quite the case and 
we find for the ratio of branching ratios with a recoiling meson X and apart from phase 
space corrections, 

R x = B ( X ^ _ tan 2 (6» 4 - 6 V ) = 4.2 x 10" 3 or 1.5 x 10~ 4 , (45) 
B(Xlo) 

for the quadratic (8 V = 39°) or linear (9 V = 36°) Gell-Mann-Okubo mass formula (Barnett, 
1996). 

The branching ratios for pn annihilation into ir~(ft and tt~uj have been measured in 
deuterium bubble chambers. The ratio R^- lies in the range 0.07 to 0.22 indicating a 
strong violation of the OZI rule (for a review, see Dover and Fishbane (1989)). The As- 
terix experiment at LEAR has measured (ft production in pp annihilation into n°(ft, rjcft, p°(ft 
and to(ft in gaseous hydrogen at NTP (50% P-wave annihilation) and in coincidence with 
atomic L X-rays (P-wave annihilation). The branching ratios for pure S-wave were then 
obtained indirectly by linear extrapolation (Reifenrother, 1991). With the corresponding 
oj branching ratios, then available from literature, the authors reported a strong violation 
of the OZI rule, especially for ir°(ft. Some of the branching ratios from Crystal Barrel are, 
however, in disagreement with previous results. We shall therefore review the direct mea- 
surement of (ft production in liquid from Crystal Barrel and then reexamine the evidence 
for OZI violation with the two-body branching ratios listed in Table [2|. 
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6.1 Annihilation into 7r°0, rjcj) and 70 

Crystal Barrel has studied the channels 

pp -» K s K l tt° and K s K L rj, (46) 

where ifs decays to 7r°7r° and r\ to 77, leading to six photons and a missing (undetected) 
Kl (Amsler, 1993d). The starting data sample consisted of 4.5 x 10 6 0-prong annihilations. 
By imposing energy and momentum conservation, the masses of the three reconstructed 
pseudoscalars and the K$ mass, a (5C) kinematic fit was applied leading to 2,834 KsK^ir 
and 72 KsK^rj events with an estimated background of 4%, respectively 36%. 

The KsKlit Dalitz plot is shown in Fig. 16. One observes the production of K*(892) 



(— * Kir) and <f){— * KsKl)- A Dalitz plot analysis was performed with the method de- 
scribed in section |7.2j . Since the C-parity of K$Kl is negative^ the contributing initial 
atomic S-state is 3 Si. One obtains a good fit to the Dalitz plot with only two amplitudes, 
one for K*K (and its charge conjugated KK* which interferes constructively) and one 
for ir°(f) with a relative contribution to the KsK^tt channel of 

K*K + KK* , , 
1 = 2.04 ±0.21. (47) 

7T U <p 

The final state K$Klt} is much simpler since only rj(fi contributes (Amsler, 1993d). 
One obtains by comparing the intensities for ir°(j) and rjcj) 

taking into account the unobserved decay modes of the r\ meson. 

Events with a Kl interacting in the Csl crystals have been removed by the selection 
procedure which required exactly six clusters in the barrel. The last two results therefore 
assume that the interaction probability for Kl in Csl does not vary significantly with 
K^ momentum. Therefore, this interaction probability needs to be determined to derive 
absolute branching ratios for n°(f) and rjcf). This number cannot be obtained directly 
by Monte Carlo simulation due to the lack of data for low energy Kl interacting with 
nuclear matter. With monoenergetic (795 MeV/c) Kl from the channel pp — > KsKl, 
Amsler (1995c) finds an interaction probability of (57 ± 3) % in the Csl barrel. This 
leads to a branching ratio of (9.0 ± 0.6) xl0~ 4 for KsKl, in agreement with bubble 
chamber data (Table |2[). 

An average interaction probability of (54 ±4)% was measured with the kinematically 
well constrained annihilation channel 2fg(— ► 7r + 7r _ )i^i7r° (Abele, 1997e). However, in 
Amsler (1995c) a somewhat lower probability was used. Updating their ir°(j) branching 
ratio one finds together with their compatible result from n ((f) — ► K + K~): 

B(pp -> ttV) = (6.5 ± 0.6) x 10" 4 . (49) 



From Eq. (48) one then obtains 

B(pp -> n<f>) = (7.8 ±2.1) x 10~ 5 . (50) 

Both numbers are slightly higher than from indirect data in gaseous hydrogen which were 
extrapolated to pure S-wave annihilation (Reifenrother, 1991), see Table ||. 



5 Note that K°K° recoiling against ir° appears as KsKs + KlKl (J pc = even ++ ) from 1 So and as KsKl 
(J PC = odd—) from 3 Si. 
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Radiative annihilation into <j> mesons should be suppressed by both the OZI rule and 
the electromagnetic coupling. Crystal Barrel has studied the channel 7^ with the reactions 



pp — > KgK^ and K + K 7, (51) 

(Amsler, 1995c). In the first reaction K$ decays to Kl is not detected and thus the 

final state consists of five photons. The KsKlJ final state was selected from 8.7 x 10 6 0- 
prong annihilations by performing a 4C fit, imposing energy and momentum conservation, 
the masses of the two pions and the K$ mass. The background reaction pp — ► KgKi 
with an interacting Kl faking the fifth photon could easily be suppressed since Kl and 
Kg are emitted back-to-back. The experimental KsKlJ Dalitz plot is dominated by 
background from KsKltt with a missing (undetected) low energy 7 from ir° decay and 



is therefore similar to the one shown in Fig. 16. The background contribution to 7^, 
mainly from ir (j) with a missing photon, was estimated by Monte Carlo simulation and 
by varying the photon detection threshold. This led to 211 db 41 7^ events corresponding 
to a branching ratio -B(7</>) = (2.0 ± 0.5) x 10 -5 , after correcting for the (updated) Kl 
interaction probability. 

A sample of 1.6 x 10 6 2-prong annihilations was used to select the second reaction in 



(51). After a cut on energy and momentum conservation (assuming kaons) the measured 
ionization loss in the JDC was used to separate kaons from pions. Again, the dominating 
background to jcp arose from 7r°0 with a missing 7. This background was subtracted by 
varying the 7 detection threshold and by keeping only 7's with energies around 661 MeV, 
as required by two-body kinematics. The signal of 29 events led to a branching ratio 
B{l<fi) = (1-9 =b 0.7) x 10~ 5 which is less precise but in good agreement with data from 
the neutral mode. The average is then 

B(pp^<fa) = (2.0 ±0.4) x HT 5 , (52) 

which updates the result from Amsler (1995c). 

6.2 (f)/u ratio 

Table |^ and Fig. 17 show the phase space corrected ratios Rx (Eq. (|45[)). The nearest 



threshold dominance factors (Eq. (|20|)) have been used but the measured ratios do not 
differ significantly from Rx- Phase space factors of the type ( |28|) lead to even larger 
ratios. For X = 7,7r°,?7 we used Crystal Barrel data. For X = to we used for ww the 
branching ratio from Crystal Barrel (multiplied by two for identical particles) and for 
uj(f) the branching ratio from Bizzarri (1971) (Table ^). For completeness we also list the 
result for X = p from Reifenrother (1991) and Bizzarri (1969) and the recent Obelix data 
for R w - (Ableev, 1995) and R w + (Ableev, 1994) in pn and np annihilation. 

Annihilation into u;7r 7r° and </>(—> KlKs)^^ can be used to extract the ratio R a 
(Spanier, 1997), where a stands for the low energy (tttt) S-wave up to 900 MeV (section 
|8.4| ). Finally, the ratio R n + n - was also measured in pp at rest (Bertin, 1996). Tableland 
Fig. [l7] show their result for 7r + 7r~ masses between 300 and 500 MeV. 

The production of 4> mesons is enhanced in all channels except r/(ft and is especially 
dramatic in 70 (from 1 5'o) and ir°(f) (from 3 Si). Several explanations for this effect have 
been proposed: 

(i) High energy reactions reveal the presence of sea ss pairs in the nucleon at high 
momentum transfers but valence ss pairs could enhance the production of <j> mesons 
already at small momentum transfers. Figure [l8| shows the OZI allowed production of ss 
mesons through the shake-out mechanism (a) and through the OZI allowed rearrangement 
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process (b) (Ellis, 1995). The fraction of ss pairs in the nucleon required to explain the 
measured ir°(j) rate lies between 1 and 19%. Deep inelastic muon scattering data indicate 
an ss polarization opposite to the spin of the nucleon (Ellis, 1995). For annihilation 
from the 3 SI state the wave function of the ( 3 Si) ss would match the wave function of 
the <f> in the rearrangement process of Fig. |l8](b), leading to an enhanced production of 
(j> mesons. The absence of enhancement in r/(fi could be due to destructive interference 
between additional graphs arising from the ss content of rj. However, this model does not 
explain the large branching ratios for the two-vector channels p°(p, uj(p and especially jcp 
which proceed from the 1 Sq atomic state. 

In the tensor nonet, the mainly ss meson is 7^(1525) and /2(1270) is the mainly uu+dd. 
Using annihilation into K^K^ir and 37r° (section ||) Crystal Barrel has measured the ratio 
of ^(1525)7r to / 2 (1270)vr° from ^o- After dividing by W (Eq. @) one finds with the 
most recent branching ratio for _/2(1525) decay to KK from Barnett (1996), see Table |l3| 
below: 

B{ VP ^m^) = _ = io _ 2 _ 

£(pp^/ 2 (1270)7r°) 

The mixing angle 9t in the 2 ++ nonet is found to be (26.1^'g) , in good agreement with 
the linear (26°) or quadratic (28°) mass formulae (Barnett, 1996). There is therefore no 
OZI violating ss enhancement from /2(1525)7r° in liquid hydrogen. 

(ii) Dover and Fishbane (1989) suggest that the 7r°0 enhancement is due to mixing 
with a four-quark state (ssqq) with mass below 2m p (Fig. |l^(c)). This exotic meson 
would then have the quantum numbers of the pp initial state ( 3 Si = 1 with 1 = 1). 
This would also explain why rjcj) (I = 0) is not enhanced and why n (j) is weak from 
1 P\ (Reifenrother, 1991). A 1 state, C(1480) — > ircf), has in fact been reported in 
ir~p — ► TT (pn (Bityukov, 1987). This isovector cannot be qq since it would decay mostly 
to TTco, which is not observed. Indeed, Crystal Barrel does not observe any 7r°u; signal 
in this mass region in pp — ► ir°ir°Lj at rest (Amsler, 1993a). Also, C(1480) has not been 
observed in pp annihilation at rest into ir + ir~(J) (Reifenrother, 1991) nor into 7r°7rV (Abele, 
1997f). In any case the large 7^ signal with the "wrong" quantum numbers 1 So remains 
unexplained by this model. 

(hi) We have seen (Eq. (17)) that the KKtt final state is dominated by K*K production 
which proceeds dominantly from the 1 = 1 3 Si state (see section |3.8j ). The (f> enhancement 
could then be due to K*K and pp rescattering (Fig.^(d)). In Gortchakov (1996) K*K 
and pp interfere constructively to produce a ir°(J) branching ratio as high as 4.6 x 10~ 4 , 
nearly in agreement with experimental data. In Locher (1994) and Markushin (1997) the 
large 7^ branching ratio simply arises from VDM: The channels p°u> and ojuj interfere 
destructively in Eq. ( j30l) (thereby lowering the branching ratio for 70;) while p°(f) and 
u>(f> interfere constructively (thus increasing the branching ratio for jcp). This conclusion, 
however, depends strongly on the phase space correction: Prescription (|^) leads to a 7^ 
branching ratio which exceeds the OZI prediction by a factor of ten (Table ||). Also, as 
pointed out by Markushin (1997), the large p°cf) an d t^<p rates remain unexplained but 
could perhaps be accommodated within a two-step mechanism similar to the one shown 
in Fig. [11(d). 

The origin of the 4> enhancement is therefore not clear. If / = 1 also dominates pp an- 
nihilation into K*K from P-states in gaseous hydrogen, then the rescattering model would 
presumably conflict with the weak ircj) production observed from 1 Pi states (Reifenrother, 
1991). With strange quarks in the nucleon, the 3 S'i contribution and hence the contri- 
bution from ss pairs in the nucleon will be diluted by the large number of partial waves 
at higher p momenta. Hence <fi production should decrease with increasing momentum. 
Also, the small r]<j) rate, possibly due to destructive interference, implies that r]'(p should 
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be abnormally large (Ellis, 1995). Finally, since /2(1525) is a spin triplet meson one would 
expect a strong production of /2(1525)7r° from triplet pp states at rest, hence 3 Pi (Ellis, 

1995) . The analysis of Crystal Barrel data in gaseous hydrogen and in flight will hopefully 
contribute to a better understanding of the <f> enhancement in hadronic reactions. 

7 Meson Spectroscopy 
7.1 Introduction 

Mesons made of light quarks u, d, s are classified within the qq nonets of SU(3)-flavor. 
The ground states (angular momentum L = 0) pseudoscalars (0 '") and vectors (1 ) 
are well established. Among the first orbital excitations (L = 1), consisting of the four 
nonets ++ , 1 ++ , 2 ++ , l + ~, only the tensor (2 ++ ) nonet is complete and unambiguous 
with the well established a 2 (1320), / 2 (1270), i^ 1525 ) and #1(1430) but another tensor, 
/2(1565) was discovered at LEAR in the 1500 MeV mass range (May, 1989). 

Before Crystal Barrel three scalar (0 ++ ) states were already well established: ao(980), 
/o(980) and ifg(1430). Further candidates have been reported and we shall discuss the 
scalars in more details below. In the 1 ++ nonet two states compete for the ss assignment, 
/i(1420) and /i(1510). In the V nonet the ss meson is not established although a 
candidate, /ii(1380), has been reported (Aston, 1988a; Abele 1997f). Many of the radial 
and higher orbital excitations are still missing. Recent experimental reviews on light quark 
mesons have been written by Blum (1996) and Landua (1996) and theoretical predictions 
for the mass spectrum can be found in Godfrey and Isgur (1985). 

Only overall color-neutral qq configurations are allowed by QCD but additional col- 
orless states are possible, among them multiquark mesons (q 2 q 2 , q 3 q 3 ) and mesons made 
of qq pairs bound by an excited gluon g, the hybrid states (Isgur and Paton, 1985; Close 
and Page, 1995). The 2 ^ state r/2(1870) has been reported by Crystal Barrel (Adomeit, 

1996) with decay rates to ^(1320)^° and / 2 (1270)77 compatible with predictions for a 
hybrid state (Close and Page, 1995). Hybrids may have exotic quantum numbers, e.g. 
1 h , which do not couple to qq. An isovector state, /3(1405), with quantum numbers 1 
has been reported (Aide, 1988a; Thompson, 1997). However, lattice QCD predicts the 
lightest hybrid, a l+~, around 2000 MeV (Lacock, 1997). 

A striking prediction of QCD is the existence of isoscalar mesons which contain only 
gluons, the glueballs (for a recent experimental review, see Spanier (1996)). They are a 
consequence of the non-abelian structure of QCD which requires that gluons couple and 
hence may bind. The models predict low-mass glueballs with quantum numbers ++ , 2 ++ 
and h (Szczepaniak, 1996). The ground state glueball, a ++ meson, is expected by 
lattice gauge theories to lie in the mass range 1500 to 1700 MeV. The mass of the pure 
gluonium state is calculated at 1550 ± 50 MeV by Bali (1993) while Sexton (1995) predicts 
a slightly higher mass of 1707 ± 64 MeV. The first excited state, a 2 ++ , is expected around 
2300 MeV (Bali, 1993). 

Since the mass spectra of qq and glueballs overlap, the latter are easily confused with 
ordinary qq states. This is presumably the reason why they have not yet been identified 
unambiguously. For pure gluonium one expects couplings of similar strengths to s's and 
uu + dd mesons since gluons are flavor-blind. In contrast, ss mesons decay mainly to 
kaons and uu + dd mesons mainly to pions. Hence decay rates to -kit, KK, rjrj and rjrj' 
can be used to distinguish glueballs from ordinary mesons. However, mixing with nearby 
qq states may modify the decay branching ratios (Amsler and Close, 1996) and obscure 
the nature of the observed state. Nevertheless, the existence of a scalar gluonium state, 
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whether pure or mixed with qq, is signalled by a third isoscalar meson in the ++ nonet. 
It is therefore essential to complete the SU(3) nonets in the 1500 - 2000 MeV region and 
to identify supernumerary states. The most pressing questions to be addressed are: 

1. What are the ground state scalar mesons, in particular is /o(980) the ss state and is 
ao (980) the isovector or are these states KK molecules (Weinstein and Isgur, 1990; 
Close, 1993) in which case the nonet members still need to be identified? Where are 
the first radial excitations and is there a supernumerary 1 = scalar in the 1500 
MeV region? Is /j(1710) scalar or tensor? 

2. In the sector, are 77(1295) and 77(1440) the two isoscalar radial excitations of rj 
and rj or is 7/(1440) a structure containing several states (Bai, 1990; Bertin, 1995), 
in particular a non-qq state around 1400 MeV? 

3. Where are the hybrid states? Is 772(1870) a hybrid and does /5(1405) really exist? 

Before reviewing the new mesons discovered by Crystal Barrel and providing clues to some 
of these issues, we shall recall the mathematical tools used to extract the mass, width, 
spin and parity of intermediate resonances in pp annihilation at rest. 



7.2 Spin-parity analysis 

The Crystal Barrel data have been analyzed using the isobar model in which the pp 
system annihilates into N "stable" particles (vr 1 * 1 , K^, K°, 7T°, r], if) through intermediate 
resonances. The decay chain is assumed to be a succession of two-body decays a — > be 
followed by b — » 6162 and c — > c\c%, etc. Final state rescattering is ignored. We shall 
calculate from the N momentum vectors the probability w o that the final state proceeds 
through a given cascade of resonances. The final state may be from real data or from 
phase space distributed Monte Carlo events to be weighted by wd- 

The spins and parities of intermediate resonances are determined using the helicity 
formalism developed by Jacob and Wick (1959) or the equivalent method of Zemach 
tensors (Zemach, 1964, 1965). Here we describe briefly the helicity formalism. Suppose 
that a mother resonance with mass thq and spin J decays into two daughters (spins Si 
and S2) with total spin S and relative angular momentum L. As quantization axis we 
choose the flight direction of the mother. The decay amplitude is given by the matrix 
(Amsler and Bizot, 1983) 

A Xl ,Xr,M = D{ M (e^)(JX\LS0X)(SX\SiS 2 Xi,-X 2 ) x BW L (m), (54) 

where the row index A = Ai — A2 runs over the (2Si + 1)(2S , 2 + 1) helicity states and the 
column index M over the 2 J + 1 magnetic substates; 6 and (ft refer to the decay angles in 
the mother rest frame. BWL{m) is the Breit-Wigner amplitude^ 

BW L (m) = — %°E?-——F L (p) y /p, (55) 



rH = ^^. (56) 



where 

m p Fl(p) 
™ Po Fl(p )' 

The mass and width of the resonance are mo and To, p is the two-body decay momentum 
and po the decay momentum for m = tjiq . The damping factors Fl (p) are given in Table 



6 The phase space factor J~p = yj2p/m should be dropped in Eq. (^5|) when events are drawn by Monte 
Carlo simulation, already assuming phase space distribution. 
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Q The matrix D is given by 



D J XM (e,<t>) = e iM 't>4 M (e) 



(57) 



where the matrix d XM {9) is the usual representation of a rotation around the y-axis, see 
for example Barnett (1996). 

We shall describe the annihilation to the observed final state by a product of matrices 
A for successive decays in the cascade. Hence we first calculate from the iV final state 
momentum vectors the angles 6 and </> for all resonances through a series of Lorentz 
boosts, apply Eq. (54) to each decay and obtain the total amplitude through matrix 
multiplications, for example 



A = [A(c -> cic 2 ) 8) A(b -> bib 2 )]A(pp -» be), 



(58) 



where ® denotes a tensor product. The matrix A has as many rows as the total final state 
spin multiplicity and has 2 J + 1 columns, where J is the total spin of the pp system. We 
define the quantization axis as the direction of one of the daughters in the first decay (the 
annihilating pp atom) for which we choose = 0,0 = and BWiCZmp) = Fl{p)\/P- 

Several decay chains of intermediate resonances may lead to the same observed final 
state of N stable particles. The transition probability W£> for chains starting from the 
same atomic state is given by the coherent sum 



WD 



w x e 



X Tr t(E a 3 A MY. «fe4)l = W X £ X Tr I E a 3 A 3 

j 



k 



(59) 



where the sums extend over all decay chains labelled by the matrices Aj. We have as- 
sumed that the initial spin-density matrix p is unity since the pp system is unpolarized. 
The phase space w and the detection probability e will be ignored for Monte Carlo events 
drawn according to phase space and submitted to the detector simulation, since w = 1 and 
e = 1 or for every Monte Carlo event. The parameters aj = aj exp(—i(f)j) are unknown 
constants to be fitted and one phase, say c/>o, is arbitrary and set to zero. For chains de- 
caying into the same resonances but with different electric charges (e.g. p + vr _ , p~vr + , p°7r° 
— > 7r + 7r~7r ) these constants are given by isospin relations. The contributions from differ- 
ent atomic states are given by incoherent sums, i.e. by summing weights wr> of the form 

©. 

As an example, let us derive the weight wrj for the annihilation channel pp — > p°p° — > 
27T + 27r - from the atomic state 1 S'o (J PC = *"). Parity, C-parity and total angular 
momentum conservation require for p°p° that L = 1 and 5 = 1 (see section 4.1). The 
first Clebsch-Gordan coefficient in (|54] ) is, apart from a constant, 

(60) 



(0A|110A)(lA|llAi,-A 2 ) =X 1 5 XlX2 , 



and hence the amplitude vanishes unless the p's are emitted with the same helicity Ai = 
A 2 7^ 0. With J = 0, the matrices (|57|) are unity and therefore A(pp) is a column-vector 
with 9 rows 

/ 1 \ 





A(pp) 



Fi{p)y/p- 



(61) 



v - 1 J 
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For p — > tt + it one finds with S = and J = 1 that L = 1 and the Clebsch-Gordan 
coefficients are unity. Hence we get with A = the 3-dimensional row-vector 

A(p) = [D 1 01 (e,<l>),D 1 00 (OA),D^ 1 (e, ( t>)}BW 1 (m). (62) 

With Eq. (|5g|) one then obtains 

A = [A( Pl ) ® A(p 2 )]A(pp) 

= [D 1 01 (d u fa)D 1 01 (d 2 , fa) - Dl_ 1 (0 1 ,<h.)Dl_ 1 (e 2 , fa)} 

xBW l (m 1 )BW l {m 2 )F 1 (p)^p, 
= isin0isin0 2 sin(0i + fa)BW 1 (m 1 )BW 1 (m 2 )F 1 (p) y /p~, (63) 

and therefore 

w D = sin 2 6>i sin 2 2 sin 2 (0i + fa)\BW 1 (m 1 )BWi [m 2 )\ 2 Ff(p)p. (64) 

The angles refer to the directions of the pions in the p rest frames, with respect to the 
flight direction of the p's. Therefore the most probable angle between the planes spanned 
by the two dipions is 90°. This angular dependence is familiar in parapositronium (0 *~) 
annihilation or 7r° decay where the 7 polarizations are preferably orthogonal (fa + fa = 
90°). However, there are two ways to combine four pions into p°p° and therefore the final 
weight is actually given by the coherent sum ( |59[ ) of two decay chains with a\ = a 2 . 

As another example of symmetrization let us consider pp annihilation into irKK which 
will be discussed in detail below. The amplitudes for annihilation through the intermediate 
K* are related through isospin Clebsch-Gordan coefficients (see for example Conforto 
(1967) or Barash (1965)). In general, annihilation may occur from 1 So or 3 Si with isospin 
I = or 1. For example, -k®K + K~ proceeds through K* + — > it°K + or K*~ — ► tt K~ 
with coefficients ol\ and a 2 equal in absolute magnitude and the two chains interfere. 
Table ||] gives the relative sign between ot\ and a 2 . Note that for 3 5*i the matrix (|57]) 
flips the sign so that the observed interference pattern is the same for 1 Sq or 3 5i, namely 
constructive in tt°K + K~, vr K°^o anc } (J = 0) ir^K^K , and destructive in (J = 1) 
ir^K^K . The signs given in Table || also apply to KK intermediate states with isospin 
i = 1 from I(pp) = and I(pp) = 1 (the latter only contributing to tt^K^K ) and for 
states with isospin i = from I(pp) = 1 (to which ir^K^K does not contribute). 

The procedure to analyze data is as follows: Phase space distributed Monte Carlo 
events are generated, tracked through the detector simulation and submitted to the re- 
construction program. As already mentioned, this procedure automatically takes care of 
the factors p, w and e. For pp annihilation into three stable particles, there are only two 
independent kinematic variables. One usually chooses the invariant masses squared m\ 2 
and 77i 2 3 . The two-dimensional distribution (Dalitz plot) is then uniformly populated for 
phase space distributed events. The procedure consists in calculating wd for each Monte 
Carlo event and to vary the constants a%, the widths, masses and spin-parity assignments 
of the resonances until a good fit to the observed Dalitz plot density is achieved. Reso- 
nances with spins larger than 2 are too heavy to be produced in pp annihilation at rest 
and are therefore ignored. 

For the Dalitz plot fits it is convenient to factorize wd in terms of the real constants 
dj and fa: 

w D = ^2 afQu + 2^2aiajRe[Qij] cos(fa - fa) + 2 ^ a,iajIm[Qij] sm(fa - fa), (65) 

i i<j i<j 

where 

Qij = Tr [AiA]). (66) 
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Dalitz plots weighted by Q„, Re[Qjj] and Im[Qij]are then produced for each pair of chains 
i,j (i < j)- The Qij, and correspondingly the weights wd, are normalized to the total 
number of real events Nt- 

Qij Qij/\ffifj, (67) 

with 

fi = J2Qn/ N T, (68) 
where the sum runs over all Monte Carlo events. One then divides the Dalitz plots into 
cells, adds them according to Eq. ( |65[ ) and builds the \ 2 

\ - (n - wu ^ 2 
^ n + w 2 D /n M c 



x A = y. , ( 69 ) 



where the sum extends over all cells. The number of real events in each cell is denoted by 
n and the number of Monte Carlo events by umc- 

With limited statistics or for more than two degrees of freedom (final states with 
more than three stable particles) the \ 2 minimization may be replaced by a likelihood 
maximization. One minimizes the quantity S = - 2 InL or 

J2 Wi [MC] J - 2 ^ Inwi [DAT] , (70) 

i=l J i=l 

where W{[MC] and Wi[DAT] are weights wry calculated for Monte Carlo and data events, 
respectively. The sums run over Nt data events and Mt Monte Carlo events. 

From the best fit the fractional contributions of the resonances in chain i are given by 

■2 



n = (71) 



where, obviously, X) r « = 1- This is a somewhat arbitrary definition which may not agree 
with the directly visible Dalitz plot densities, because interferences beween the chains are 
neglected in Eq. (|7T|) . One may define alternatively 



ri = c4, (72) 

but then Yl r i may differ significantly from unity in the presence of strong interferences. 
Hence decay branching fractions for broad interfering resonances are not measurable un- 
ambiguously. This is an unavoidable caveat to keep in mind when extracting the internal 
structure of broad states from their decay branching ratios. 



7.3 K- matrix analysis 

The Breit-Wigner factors ( |55| ) violate unitarity when two resonances with the same quan- 
tum numbers overlap and decay into the same final state. Also, they do not describe 
distortions in the mass spectrum that occur around kinematical thresholds. For example, 
the /o(980) — * tttt appears as a dip rather than a peak in the irn mass spectrum of elastic 
■kit scattering, due to the opening of the decay channel /o(980) — ► KK (Au, 1987). 

This behaviour can be described with the K-matrix formalism. A detailed description 
can be found in Chung (1995) and I shall only recall the formulae used in the analysis of 
Crystal Barrel data. Consider for instance the four scattering reactions 

■KIT — > KK KK — > KK 

KK — > ttk KK^KK 
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The transition amplitude T for a given partial wave is described by the 2x2 ivT-matrix 

T = (1 — iKp)~ x K (74) 
with the real and symmetric matrix 

Kij(m) = ° m 2 3 _ ° 2 ° ^( m )^ ( m ) + ( 75 ) 

The sum runs over all resonances with iC-matrix poles m a decaying to tttt and KK with 
(real) coupling constants 7 ai and 7 Q2 , respectively, where 

7^+7^ = 1- (76) 
The factors Z? ai are ratios of Blatt-Weisskopf damping factors (Table ||) 

where L is the decay angular momentum, p the ir or K momenta and p ai their momenta at 
the pole mass m a . The optional real constants Cj,- allow for a background (non-resonating) 
amplitude^. In Eq. ( |74| ) the matrix p(m) describes the two-body phase space and is 
diagonal with p\\ = p\ = 2p- K /m and P22 = P2 = 2pi</ m - For masses far above kinematical 
threshold pi ~ 1 and below iTX threshold /?2 becomes imaginary. 

The (if-matrix) partial width of resonance a to decay into channel i is defined as 

T ai (m a ) =-flT' a pi(m a ), (78) 

and the (iT-matrix) total width as 

r a = £r ai . ( 79 ) 

i 

For a resonance with mass far above kinematic thresholds one obtains the partial and 
total widths 

r Qi = 7 ^r' a , r Q = r' a) (so) 

respectively. For one resonance and one channel (elastic scattering) the i^-matrix reduces 
to 

K = (8!) 

and T (Eq. ([74])) reduces to the relativistic Breit-Wigner 

rr _ m ToB 2 (m)/p(mo) 

1 — — 2 2 : 77 — V' v 82 ** 

rrig — m z — imr)l (m) 

with 

r(m) = ro-^UV)- (83) 

For a resonance far above threshold and with To << mo we get the familiar expression 

r /2 



mo — m — iro/2 



(84) 



7 For the tttt S-wave a factor (to 2 — 2ni^)/m 2 is multiplied to the if -matrix to ensure a smooth behaviour 
near threshold. 
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Normally, the mass m R and width Tr of a resonance are obtained from the poles of the 
T-matrix. Extending the mass m to complex values we find from (|84|) the poles at 



m P = m R -i^Y, (85) 

with rriR = itiq and Tr = Tq. In general, however, uir does not coincide with the pole 
of the .fT-matrix and Tr is different from the K-matrix width. For example, for two non- 
overlapping resonances far above threshold, the .fT-matrix, T-matrix and Breit-Wigner 
poles coincide. As the resonance tails begin to overlap the two T-matrix poles move 
towards one another (for an example, see Chung, (1995)). 

The -kit S-wave scattering amplitude is related to the -kit phase shift 5 and inelasticity 
rj through the relation 

pi(m)Tii(m) = . (86) 

li 

According to Eq. (^) the corresponding X-matrix then reads for pure elastic tttt- 
scattering (r/ = 1) 

y r , s tan 5(m) 

Ku m = (87) 

pi{m) 

and becomes infinite at m = mo, when 5 passes through 90°. However, the amplitude 
T does not, in general, reach a resonance when 5 = tt/2. As an example, consider the 
tttt S-wave scattering amplitude described by the amplitude (|86|) in the complex plane 
(Argand diagram): The intensity |T| 2 reaches its maximum value around 850 MeV (5 = 



90°), loops back and passes rapidly through the KK threshold (see Fig. 23 below and 
Au (1987)). At ~1000 MeV |T| 2 reaches its minimum value (5 = 180°) and then starts a 
new (inelastic) loop. The /o(980) then appears as a hole in the tttt intensity distribution. 
We shall return to the S-wave Argand diagram when discussing the fits to Crystal Barrel 
data. 

Consider now the production of a resonance a in pp annihilation. In the isobar model, 
the resonance is assumed not to interact with the recoiling system. The coupling strength 
to pp is denoted by the (complex) constant (3 a while 7 Qi describes its decay strength into 
channel i (say tttt for i = 1 and KK for i = 2). Following Aitchison (1972) the amplitudes 
are given by the components of the vector 

T = (1 - iKp)~ 1 P. (88) 

The isT-matrix now describes the propagation of the channel i through the resonances a 
while the vector P describes their production. P and K share the common poles m a so 
that T remains finite at the poles. The vector P is given by 

P 3 ( m ) = 2 1 ' ( 89 ) 

a m a ~ m 

where the sum runs over all resonances. For a single resonance feeding only one decay 
channel we again obtain from Eq. ( ]88| ) a Breit-Wigner distribution of the form (^) with 
coupling strength (3: 

j3m T B(m)l >(m ) 
m\ — m 2 — imoT(m) 

Let us now assume a series of resonances with the same quantum numbers decaying 
into two final states. The amplitude for the first final state is given by Eq. (]88|): 



^ _ (1 - iK 2 2P2)P\ + 1K12P2P2 / gl x 
1 - P1P2D - i(piK u + £2-^22) ' 
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with 

D = K U K 22 - Kf 2 . (92) 
As an example, consider a single resonance, say ao(980) decaying to r/n and KK. In this 



case D = and B{m) = 1 (S-wave). We then obtain from Eq. (91) the formula (Flatte, 
1976) 

T(mr) = — = — b91 (93) 



m 2 _ m 2 _ i{pig 2 + ml) 



with 

9i = It 

and 

The phase space factors are 



m T' Q J2g! = m T' , (94) 



b = (3,m T' . (95) 



2 2 / 4m 2 
pi(m) = — prj and j0 2 (m) = — pk = A/1 (96) 

TO TO V TO/ 



The T(KK) amplitude is obtained by interchanging the labels 1 and 2 in Eq. fl93|) . Below 
XX threshold /92 becomes imaginary. Compared to pure r]ir decay this leads to a shift of 
the resonance peak and to a narrower and asymmetric distribution of the observed signal 



in the r]ir channel. This is shown in Fig. |19| for g\ = 0.324 GeV, g 2 = 0.279 GeV (hence 



T'q = 0.43 GeV). These parameters have been extracted from the ao(980) contribution 
to pp — > i]tt°tt° and KKir (section |8.8[) . A width of 54.12 ± 0.36 MeV was determined 
directly from the ao(980) — ► t]tt signal in the annihilation channel ior]ir° (Amsler, 1994c), 



in good agreement with the observed width in Fig. |19[ Also shown in Fig. |19| is the 
expected distribution for T' Q = 0.43 GeV, assuming no KK decay. The observed width 
Tq = r' pi(mo) increases to 0.28 GeV. 

The standard procedure in the analysis of Crystal Barrel data is to replace the Breit- 



Wigner function (|55|) by T (Eq. (|8q) ) and to fit the parameters g ai , (3 a and m a . The 
resonance parameters rriji and Tr are then extracted by searching for the complex poles 
(Eq. (p5|)) of the matrix T. A one-channel resonance appears as a pole in the second 
Riemann-sheet and a two-channel resonance manifests itself as a pole in the second or 
third Riemann-sheet (Badalyan, 1982). 

Some of the Crystal Barrel Dalitz plots have also been analyzed using the N/D for- 
malism (Chew and Mandelstam, 1960) which takes into account the direct production 
of three mesons and also the final state interaction. In this formalism the amplitude T 
has the same denominator as, e.g., in Eq. (pi]), but the numerator allows for additional 
degrees of freedom (Bugg, 1994). 



8 Annihilation at Rest into Three Pseudoscalars 

Proton-antiproton annihilation at rest into three pseudoscalars is the simplest process to 
search for scalar resonances + — > 0~0~, the recoiling third pseudoscalar removing the 
excess energy. The annihilation rates for these processes in liquid hydrogen ( 1 5o atomic 
state) are reasonably large since no angular momentum barrier is involved. Channels with 
three pseudoscalars have been studied earlier, but essentially in the 2-prong configuration 
and with limited statistics, e.g 7r + 7r~7r° with 3,838 events (Foster, 1968b), ir + ir^r] with 
459 events (Espigat, 1972) and tt + tt~i]' with 104 events (Foster, 1968a). The samples 
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collected by the Asterix collaboration at LEAR (May, 1989; Weidenauer, 1990) are larger 
but were collected from atomic P-states. The ir + ir~ir° final state revealed the existence of 
a tensor meson, /2(1565), produced from P-states (May, 1989, 1990). In the kaonic sector, 
data were collected in bubble chambers for the final states tt^K^Ks (2,851 events) and 
ir°KsKs (546 events) in the experiments of Conforto (1967) and Barash (1965). Branching 
ratios for annihilation into three mesons are listed in Table [| 

Annihilation with charged pions is dominated by p(770) production which complicates 
the spin-parity analysis of underlying scalar resonances in the tttt S-wave. Also, both 1 So 
and 3 Si atomic states contribute. All-neutral (0-prong) channels are therefore simpler to 
analyze but more complex to select due to the large 7-multiplicity. The channel 7r 7r°7r° 
with 2,100 events has been reconstructed earlier with optical spark chambers (Devons, 
1973). The existence of a scalar resonance decaying to tttt with mass 1527 and width 101 
MeV was suggested in the 37r° channel and in its ir~ir~ir + counterpart in pn annihilation 
in deuterium (Gray, 1983). This was actually the first sighting of /o(1500) which will be 
discussed below. 

The sizes of the data samples have been vastly increased by Crystal Barrel. We shall 
first review annihilation into three neutral non-strange mesons. We start from 67 final 
states and select the channels pp — > 7r°r/r/, 37r°, tt^tt !], ir°7]r]' and by assuming total 

energy and total momentum conservation and constraining the 27 masses to vr°,r/ and 
rj' decays (7C fits), excluding any other possible configuration: Events are accepted if 
the kinematic fit satisfies the assumed three-pseudoscalar hypothesis with a confidence 
level typically larger than 10%. Background from the other 67 channels is suppressed 
by rejecting those events that also satisfy any other 67 final state hypothesis (including 
the strong loco, uj — > 7r°7) with a confidence level of at least 1%. The absolute branching 
ratios for the 67 channels are determined by normalyzing on the branching ratio for 
pp —* ujuj. These three-pseudoscalar channels have all been analyzed and we now review 
the salient features in their Dalitz plots. Results on kaonic channels are appended to the 
next sections. 



8.1 pp — > 7T°TJT] 

The first evidence for two 1 = scalars in the 1400 MeV mass region, now called /o(1370) 
and /o(1500), was obtained from a reduced sample of 2.3 x 10 4 ir°r]r] events (Amsler, 
1992c). The invariant mass distributions are shown in Fig. ^0[ The two scalars decaying 
to rjrj are also observed when one rj decays to 3ir° (IO7 final state), a channel with entirely 



different systematics (Fig. 2C(b)). The distributions in Fig. po|(a) and (b) are nearly 
identical. 

An amplitude analysis of the Dalitz plot distribution for the 67 final state was per- 
formed with the method outlined in the previous section. However, Breit-Wigner functions 
of the form ( |55|) were used to describe the resonances. The fit required J PC = ++ for 
both rjr) resonances. The (Breit-Wigner) masses and widths were m =1430, V = 250 and 
m = 1560, T = 245 MeV, respectively. Note that the width of the upper state appears 
smaller in Fig. |20|(a,b), due to interference effects. This state may be identical to the 
scalar meson observed earlier by the GAMS collaboration at 1590 MeV in the rjrj and rjrj' 
mass spectrum of high energy 7rA-interactions (Aide, 1988b; Binon, 1983, 1984). 

The final analysis of this channel was performed with a tenfold increase in statistics, 
namely 3.1 x 10 4 ir°r]r] events from 0-prong data and 1.67 x 10 5 ir°r]r] events from a triggered 
data sample requiring online one ir° and two tj mesons (Amsler, 1995e). The Dalitz plot 
is shown in Fig. |2l](a). The horizontal and vertical bands are due to ao(980) decaying to 
r/7r°. One also observes diagonal bands which correspond to the two states decaying to rjrj. 
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A residual incoherent flat background of 5%, mainly due to — ► 77 with a missing 
photon, has been subtracted from the Dalitz plot before applying the amplitude analysis, 
this time with the full i^T-matrix formalism. 

Since S-wave dominates in liquid, the channel ir°r]r] proceeds mainly through the 1 iS'o 
atomic state. The ijtt S-wave was parametrized by a 2x2 i^-matrix with poles from 
ao(980) and ao(1450). The parameters were taken from the 7r°7r°r/ analysis (section |8.2| ), 
leaving the production constants (3 free. A contribution from 02 (1320) (r]7r D-wave) with 
fixed mass and width was also offered to the fit. The r\r\ S- and D-waves were described 
by one-channel X-matrices. Annihilation from atomic P-states was not included in the 
fit except for tensor mesons (e.g. the expected /2(1565)) decaying to r]rj. In fact the fit 
demands a contribution from a tensor meson with mass ~1494 and width ~155 MeV, 
produced mainly from P-states. 

The best fit was obtained with two poles for the 7777 S-wave. The resonance parameters 
(T-matrix poles) are: 

/ (1370):m = 1360 ± 35, T = 300 - 600 MeV, 

/o(1500):m = 1505 ± 15, T = 120 ± 30 MeV. (97) 

The ET-matrix mass and width of /o(1500) are 1569 and 191 MeV, respectively, in accord 
with the Breit-Wigner parameters of the GAMS resonance (Binon, 1983). This state, pre- 
viously called /o(1590), and /o(1500) are therefore assumed to be identical. The branching 
ratio for ifir\r\ is given in Table |9] and the products of resonance production and decay 



branching ratios are listed in Table 1C 



8.2 pp — > 71 ' 71 77 

This channel is relevant to search for isovector ++ states decaying to rjir. The 7r°7r°r/ 
Dalitz plot (2.8 x 10 5 events) is shown in Fig. |2~l|(b). Qualitatively, one observes do (980) 
and a2(1320) decaying to r/n and /o(980) decaying to hit. The strong interference patterns 
point to coherent contributions from a single pp atomic state ( So). 

An amplitude analysis based on the ET-matrix formalism (and, alternatively, the N/D 
formalism) has been performed, assuming pure S-wave annihilation (Amsler, 1994b). The 
-kit S-wave was described by two poles, one for /q(980), coupling to irir and KK, and one 
for /o(1370). Elastic 7T7r-scattering data (Grayer, 1974; Rosselet, 1977) were included in 
the fit. The -kit D-wave (/2(1270)) was also introduced but turned out to be negligible. 
The rjn D-wave was described by one pole for 02 (1320). The rjir S-wave was described 
by a one-pole 2x2 X-matrix for ao(980) with couplings g\ to rjir and 52 to KK. Since 
decay to KK was not measured, 52 was obtained indirectly from the r]n line shape. The 
fit yielded gi = 0.353 GeV and 

— ~ 0.88. (98) 
9i 

These amplitudes were, however, not sufficient to describe the data. A satisfactory fit 
was obtained by adding (i) a second pole to the r\n S-wave, (ii) a second pole to the rjir 
D-wave and (iii) an tjtt P-wave. Branching ratios are given in Table |l^. The branching 
ratio for 7r°7r°r/ (Table ^) is in excellent agreement with the one derived from the channel 
7r°7r°77 -> 7r + 7r"3vr (Amsler, 1994d). 

The main result was the observation of a new isovector scalar resonance in the tjtt 
S-wave: 

a (1450) : m = 1450 ±40, T = 270 ± 40 MeV. (99) 

This resonance manifests itself as a depletion in the bottom right (or top left) corner of 
the Dalitz plot (Fig. ^ (b)). Evidence for the ao(1450) decaying to t]tt is also reported 
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in an analysis of the channel 7r + 7r~ry which requires, in addition, an amplitude for prj 
production from 3 Si (Abele, 1997a). 

The rjTr D-wave contribution corresponds to a 2 ++ resonance around 1650 MeV, called 
4(1650) in Table [0J with a width of about 200 MeV. This state could be the radial 
excitation of a2 (1320). The mass and width of the structure in the r]ir P-wave (exotic 
1"+) are not well defined. They vary from 1200 to 1600 MeV and from 400 to 1000 MeV, 
respectively, without significant changes in the x 2 - The 1 *~ p(1405) reported by Aide 
(1988a) is therefore not confirmed by the present data. 



8.3 pp -> 7r 7r°7r° 

The first analysis of this channel used a sample of only 5.5 x 10 4 events and reported an 
isoscalar 2 ++ meson at 1515 ± 10 MeV with width 120 ± 10 MeV, decaying to 7r°7r° (Aker, 
1991). This state was identified with /2(1565) that had been observed before in the final 
state 7r + 7r~7r° in hydrogen gas (May, 1989, 1990). P-wave annihilation from 3 Pi and 3 i-2 
was therefore allowed when fitting the 3ir° channel. Resonances in the tttt S-wave were 
described by the tttt elastic scattering amplitude, replacing the Breit- Wigner amplitude 
by 

BW (m) = ™ (ri(m)eMmm)]-l\ 
P V 2z J 

according to Eq. d83), where 5 and rj were taken from the Argand diagram of Au (1987). 
This is an approximation assuming equal production strengths for all resonances in the 
tttt S-wave, which is reasonable for the 3ir° channel, as I will show below. 

A statistical sample an order of magnitude larger then revealed a new feature (Amsler, 
1994e) which was consolidated by a reanalysis of the early Crystal Barrel data (Anisovich, 
1994): the presence in the Dalitz plot of a narrow homogeneously populated band from 
a scalar resonance, /o(1500), decaying to 2vr°. The 3vr° Dalitz plot is shown in Fig. |2l](c) 



and the 27r° mass projection in Fig. I22. 



Qualitatively, one observes the following features: the population along the tttt mass 
band marked /2(1270) increases at the edges of the Dalitz plot indicating that one decay 
7T° iS preferably emitted along the flight direction of the resonance. This is typical of 
a spin 2 resonance decaying with the angular distribution (3cos 2 # — l) 2 from 1 So or 
(1 + 3cos 2 <9) from 3 P X . The blobs labelled / 2 (1565) at the corners correspond to an 
angular distribution sin 2 # from another spin 2 resonance produced from 3 P 2 > together 
with constructive interference from the two tttt S-waves. The /o(980) appears as a narrow 
dip in the tttt S-wave. The new feature is the homogeneous narrow band marked /o(1500) 
which must be due to a spin state. 

The analysis of the full data sample was performed with the A-matrix formalism 
(Amsler, 1995f). A 2 x 2 A-matrix with three poles was sufficient to describe the tttt 
S-wave. Elastic tttt scattering data up to 1200 MeV from Grayer (1974) and Rosselet 
(1977) were also included in the fit. The contributing scalar resonances are /o(980) and 

/ (1370) : m ~ 1330, V ~ 760 MeV, 

/o(1500) : m = 1500 ±15, T = 120 ± 25 MeV. (101) 

A 4-pole A-matrix helps to constrain the /o(1370) parameters, giving 

/ (1370) : m = 1330 ±50, T = 300 ± 80 MeV, (102) 

and the 4th pole corresponds to a 600 MeV broad structure around 1100 MeV (called 
/o(400 — 1200) by the Particle Data Group (Barnett, 1996)) and also reported by a recent 
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reanalysis of tttt S-wave data (Morgan and Pennington, 1993). The data are therefore 
compatible with 3 or 4 poles and it is not obvious that /o(1370) and /o(400 — 1200) are 
not part of the one and the same object. 

A one-dimensional /T-matrix describes the tttt D-wave. In addition to /2(1270) one 
finds 

72(1565) : m ~ 1530, V ~ 135 MeV. (103) 

The fractional contribution of P-waves is 46%. Without P-waves the fit deteriorates 
markedly but the /o(1370) and /o(1500) parameters remain stable. We shall return to 



the /2(1565) in the discussion below. Branching ratios are given in Table 10 



8.4 Coupled channel analysis 

A simultaneous fit was performed to the channels 7T°r]r], tt 7t°7] and 37r° (Amsler, 1995g) 
using the full data samples presented in the previous sections with, in addition, the tttt- 
scattering data up to 1200 MeV. However, pure S-wave annihilation was assumed. A 3 
x 3 i'T-matrix with 4 poles was used to describe the tttt S-wave coupling to tttt, rjrj and 
the at that time unmeasured KK through the resonances /o(980), /o(1370), /o(1500), 
taking also into account the broad structure around 1100 MeV. Common (5 a parameters 
(Eq. fl89|)) were introduced to describe the production of resonances associated with the 
same recoiling particle. For example, /o(1500) recoiling against 7T° is produced with the 
same strength in ir rjr] and 3n°. The tjtt S-wave was described by a 2 x 2 /Y-matrix for 
ao(980) and ao(1450). The tttt, r/r/ and r/n D-waves were treated with one dimensional 
X-matrices, the latter including a2(1320) and a 2 (1650) with pole parameters taken from 
the 7r°7r°77 analysis of section |8.2| . 

The branching ratios are given in Table [l(]. Note that the tttt S-wave includes contri- 
butions from / (400- 1200), / (980) and / (1370) (but excluding / (1500)) which cannot 
be disentangled due to interferences. The fit is in good agreement with the single channel 
analyses and constrains the resonance parameters. Hence a consistent description of all 
three annihilation channels was achieved with the following main features: 

1. The data demand three scalar resonances in the 1300 - 1600 MeV region: 

a (1450) : m = 1470 ± 25, V = 265 ± 30 MeV, 
/ (1370) : m = 1390 ±30, V = 380 ± 80 MeV, 
/o(1500) : m = 1500 ±10, T = 154 ± 30 MeV. (104) 



The broad scalar structure around 1100 MeV (/o(400 — 1200)) has very different 
pole positions in sheets II and III, making a resonance interpretation of this state 
difficult. 

The production data demand a larger width (~ 100 MeV) for /o(980) than the tttt 
scattering data alone (~ 50 MeV, according to Morgan and Pennington (1993)). 

A tensor meson is observed in the tttt D-wave with mass 1552 and width 142 MeV, 
in accord with May (1889, 1990) notwithstanding the absence of atomic P-waves 
in the present analysis. A structure is also required in this mass range in the 7777 
D-wave. This state is not compatible with 7^(1525) which is produced with a much 
lower rate, as we shall see in section ^7|. 



The 7T7T scattering amplitude T (Eq. (pa)) obtained from the fit to the elastic scattering 



data and the Crystal Barrel data is shown in Fig. 23. Note that Crystal Barrel KK data 



are not yet included and therefore the third if-matrix channel also contains by default all 
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other unmeasured inelasticities. Figure shows the 1 = S-wave production intensity 
|T| 2 (Eq. (HH)) for the three annihilation channels. The maxima around 1300 and 1550 
MeV correspond to the if-matrix poles for /q(1370) and /o(1500). It is instructive to 
compare the dip around 1000 MeV in the tttt S-wave for the 3tt° channel to the peak in 
the 7r°7r°7/ channel, both due to /o(980). This is produced by interferences between the 
tttt S-waves in 3tt° which have the opposite sign to the interference between the tttt and 
the r\TT S-waves in rjTT°TT . The tttt S-wave in tttt scattering exhibits qualitatively the same 
behaviour as in pp annihilation into 3tt°, namely sharp minima around 1000 and 1450 
MeV. The ansatz (100) used in several Crystal Barrel analyses (e.g. in Aker (1991)) for 



the tttt S-wave is therefore a rough but reasonable approximation. 

8.5 pp — > 7T 777/ 

Another piece of evidence for /o(1500) stems from TT°rjrj' (Amsler, 1994f). This channel 
was also reconstructed from the 67 final state. Since rj decays to 77 with a branching ratio 
of only 2.1% the data sample is small (977 events). Most of these events were collected 
with the online trigger requiring one 7r°, one 77 and one rj . The TT®r\rj Dalitz plot shows an 



accumulation of events at small r\rj masses. Figure 25 shows the r\rj mass projection. The 



r\rj mass spectrum from the same annihilation channel, but with rj — ► r/7r + 7r , has entirely 



different systematics but is identical (inset in Fig. 25). The enhancement at low masses is 
due to a scalar resonance since the angular distribution in the rjri' system is isotropic. A 
maximum likelihood fit was performed to the 67 channel using a (damped) Breit-Wigner 
according to Eqs. fl5"5|), ( |56| ) and a flat incoherent background. The resonance parameters 
are 

/o(1500) : m = 1545 ±25, T = 100 ± 40 MeV. (105) 

The branching ratio is given in Table [l(| The /o(1500) mass is somewhat larger than 
for 3tt° and TT rir). However, Abele (1996a) points out that a constant width in the 
denominator of the Breit-Wigner function yields a mass around 1500 MeV. This is a more 
realistic procedure since the total width at the r\rj threshold remains finite due to the 
channels tttt and r\r\. However, this does not modify the branching ratio for /o(1500) — ► r\rj . 
We shall therefore ignore Eq. ( |105| ), when averaging below the /o(1500) mass and width. 

8.6 pp — > 7r°7r°?7 / 

With 0-prong data this final state is accessible through rj — > 27 or rj — * r/7r°7r°, leading 
to 6, respectively 10 photons (Abele, 1997g). The branching ratios from both final state 
configurations agree (Table ^). A sample of 8,230 IO7 events were kinematically fitted to 
4ttV The TT°TT°r] mass distribution (Fig. ^6|) shows a sharp signal from rj and a shoulder 



around 1400 MeV due to the E meson that will be discussed in section 11. 

The TT TT rj' Dalitz plot was obtained by selecting events in the rj peak and subtracting 
background Dalitz plots from either sides of the peak. It shows a broad accumulation 
of events in its center which can be described by a dominating tttt S-wave and small 



contributions from a2(1320) and ao(1450) with branching ratios given in Table 10. These 
resonances are included in the fit with fixed mass and width. The ratios of rates for 
02(1320) and oq(1450) decays into t\tt and rjir can be predicted from SU(3) and compared 



with measurements. This will be discussed in section 10. 1] . The analysis of the 67 Dalitz 



plot (3,559 events) leads to similar results. Figure ^7] shows that ao(1450) is required for 
a satisfactory description of this annihilation channel. 
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8.7 pp tt°K l K l 

The isoscalar /o(1500) has been observed to decay into 7r°7r°, r\r\ and 777/ . To clarify its 
internal structure it was essential to also search for its KK decay mode. In a previous 
bubble chamber experiment (Gray 1983) no KK signal had been observed in the 1500 
MeV region in pp annihilation into KKtt, leading to the conclusion that the /o(1500) 
coupling to KK is suppressed (Amsler and Close, 1996). However, no partial wave analysis 
was performed due to limited statistics. 

Crystal Barrel has therefore searched for /o(1500) in the annihilation channel 7t°KlKl 
(Abele, 1996b; Dombrowski, 1996) which proceeds from the 1 Sq atomic state. All- 
neutral events were selected with three energy clusters in the barrel and the channel 
pp — ► 7t°KlKl could be reconstructed by measuring the vr°(— > 27) momentum and the 
direction of one of the Kl which interacts hadronically in the Csl barrel. The main 
background contribution was due to events with a reconstructed 7r° and one additional 
7 which happen to fulfil the tt°KlKl kinematics but for which one or more 7's have 
escaped detection. This background (~ 17%) could be removed by subtracting a Dalitz 
plot constructed from 27 pairs with invariant masses just below or above the ir° mass. 

Further background contributions were due to pp — * ojKlKl where u decays to 7r° 7 
and both Kl are not detected, and pp — ► K$Kl where one photon from K$ — ► vr°7r° — > 47 
and the Kl are undetected. These events can be removed with appropriate mass cuts. 
Background from final states like wq, wk° and 37r° were studied by Monte Carlo simulation. 
The total residual background contamination was (3.4 ± 0.5)%. 

The background subtracted Dalitz plot is shown in Fig. [H](d). This plot has not been 
symmetrized with respect to the diagonal axis since one Kl is detected (vertical axis) 
whilst the other is missing (horizontal axis). The interaction probability as a function 
of Kl momentum was measured by comparing the intensities along the two K* bands, 
K* — ► Kltt . The interaction probability was found to be flat with K l momentum, but 
increasing slowly below 300 MeV/c (Dombrowski, 1996). The tt°KlKl Dalitz plot shown 
in Fig. ^(d) is already corrected for detection efficiency and for Kl decay between the 
production vertex and the crystals. Its symmetry along the diagonal axis is nearly perfect, 
indicating that acceptance and backgrounds have been taken into account properly. 

One observes signals from if* (892) decaying to Kir and ci2(1320)//2(1270) decaying 
to KK. The accumulation of events at the edge for small KK masses is due to the tensor 
^(1525) which is observed for the first time in pp annihilation at rest. These resonances 
and the broad .Kg (1430) were introduced in a first attempt to fit the Dalitz plotQ. The 
if -matrix for the Ktt (i = 1/2) S-wave was written as 

K = -^- 2 + ^Eo/Po, (106) 
2 + abp z m,Q — m z 

where the first term describes the low energy Kit scattering (a is the scattering length 
and b the effective range) and the second term describes the if(*(1430) resonance. The 
parameters mo, To, a and b were determined by fitting the scattering amplitude T (Eq. 
@) to the Kit phase shifts (Aston, 1988b). The fit is shown in Fig. ||. Note that 
resonance occurs at 8 ~ 120°. The corresponding mass and width (T-matrix pole) for 
.K"o(1430) are m = (1423 ± 10) MeV and (277 ± 17) MeV, in close agreement with Aston 
(1988b) who used a different parametrization and found m,Q = 1429 ±6 and Tq = 287 ±23 
MeV. 



8 The absence of threshold enhancement from ao(980) or /o(980) at the upper right border of the Dalitz plot 
could be due to destructive interference between these states and/or to the loss of acceptance close to the KK 
threshold. 



35 



A one-pole .fT-matrix for a scalar resonance was added for the peak around 1450 MeV 
in the KK mass distribution (Fig. 29). The fit clearly fails to describe the KK mass 



spectrum (dashed line in Fig. |29| ), The second attempt assumed a if-matrix with two 
scalar resonances, /o(1370) and /o(1500), in the KK amplitude. The fit now provided a 
satisfactory description of the Dalitz plot and the KK mass spectrum (full line in Fig. 
|29| ). However, the isovector ao(1450) is also expected to decay to KlKl and one cannot 
distinguish between isovectors and isoscalars from the tt°KlKl data alone. Therefore a 
Breit-Wigner was added for the isovector ao(1450) with fixed mass and width taken from 
the coupled channel analysis (section |8.4j ). The /o(1370) and /o(1500) poles are stable, 
nearly independent of ao(1450) contribution. One finds 

/o(1370) : m = 1380 ± 40, T = 360 ± 50 MeV, 

/o(1500) : m = 1515 ± 20, T = 105 ± 15 MeV, (107) 

in agreement with the resonance parameters measured in the other decay channels. 

Due to uncertainties in the Kl interaction probability the branching ratio for n KlKl 
could not be determined accurately. The branching ratios given in Table [XT] (Abele, 1996b; 
Dombrowski, 1996) are therefore normalized to the known branching ratio for ir°KsKs 
(Armenteros, 1965; Barash, 1965)^]. 

The branching ratios for pp — > /o(1370) and /q(1500) — > KlKl also vary with ao(1450) 



contribution (Fig. 30). In Table |TT] the branching ratios for /o(1370) and /o(1500) are 
therefore derived from Fig. 3C assuming an ao(1450) contribution ro to the tt°KlKl 



final state, derived from its measured contribution r = (10.8 ± 2.0)% to ir K^Kl (next 
section): 

r B(pp — ► tt K^Ke) , Nrw . 

= (M*"»*. (108) 

where we have used the branching ratios from Armenteros (1965) and Barash (1965). 

8.8 pp -> -k^K^Kl 

This channel selects only isospin 1 resonances decaying to KK and therefore permits 
a direct measurement of the contribution from isovectors to KKir, in particular from 
ao(1450). Crystal Barrel has studied the reaction pp — » it K^Kl with a non-interacting 
Kl (Heinzelmann, 1996; Abele 1997e). This channel is selected from 7.7xl0 6 triggered 
2-prong data by requiring exactly two clusters in the barrel from 7r ± and K T . Particle 



identification is achieved by measuring the ionisation in the drift chamber (Fig. 31) 
and a (1C) kinematic fit ensures momentum and energy conservation. The background 
contribution, mainly from 7r + 7r~7r°, is about 2%. The Dalitz plot (Fig. 32) has been 



corrected for background and acceptance, and for the Kl interaction probability. The 
latter was determined by reconstructing the channel 7r°i^5'(^ tt + tt^)Kl with and without 
missing Kl- The branching ratio 

B(pp -> t^K^Kl) = (2.91 ± 0.34) x 10~ 3 (109) 

is in excellent agreement with the one given in Table [ll| for tt^K^Ks from bubble chamber 
experiments (Armenteros, 1965; Barash, 1965). 

The Dalitz plot shows clear signals from K*(892), a 2 (1320) and a (980). The iT(892) 
and 02(1320) were parametrized by Breit-Wigner functions. The ao(980) was described 



9 The isospin contributions from Sq to the K*K system cannot be determined in this channel since the 
bands interfere constructively for both I — and 1=1. 



36 



by a 2 x 2 i-T-matrix (Flatte formula @). The Kit S-wave (^(1430)) was treated 
using the data from Aston (1988b), as described in the previous section. In contrast to 
k^KlKl, both atomic states 1 Sq and 3 Si contribute. The 1 = and 1 = 1 contributions 
to K*(892)K can be determined from the (destructive) interference pattern around the 
crossing point of the K* bands in Fig. |3^. 

The fit, however, did not provide a satisfactory description of the Dalitz plot and the 
fitted 02(1320) mass, 1290 MeV, was significantly lower than the table value, a problem 
that had been noticed earlier in bubble chamber data (Conforto, 1967). A substantial 
improvement in the x 2 (Fig- |33| ) was obtained when introducing the ao(1450) as a second 
pole in the K-matrix, together with ao(980), leading to the resonance parameters (T- 
matrix pole) 

a (1450) : m = 1480 ± 30, r = 265 ± 15 MeV, (110) 

in agreement with the rjn decay mode. The a<i (1320) mass now became compatible with 
the table value (Barnett, 1996). We shall show below that the contribution to n K^Kl 
of (10.8 ± 2.0) % (Fig. |33|) is consistent with predictions from the 7r°7r°?7 channel, using 
SU(3). _ 

An even better fit was achieved by adding a broad structure in the KK P-wave, 
presumably from the radial excitations p(1450) and p(1700), which could not, however, 
be disentangled by the fit. 

For ao (980) one finds for the T-matrix pole in the second Riemann sheet the mass 987 
db 3 MeV and the width 86 db 7 MeV. Using the tjtt decay branching ratio from the 7r°7r°7/ 
analysis (Table |i~0| ) one also obtains the ratio 

S(a (980) -► r/vr) y J 

With the coupling gi = 0.324 GeV to ijtt (Amsler, 1994b) the coupling 02 to KK can 



be tuned to satisfy Eq. (Ill) by integrating the mass distributions over the ao(980) 
distribution (Fig. |l9|). One obtains #2 = 0.279 MeV and the ratio §2/91 = 0.86, is in 
excellent agreement with the estimate from the line shape in the 7r°7T°r] channel (Eq. (p8|)). 



The branching ratios are given in Table |11. The intermediate K*K is largest in the 
(1=1) 3 Si channel, a feature that was noticed before (Barash, 1965; Conforto, 1967) and 



that we have used in section 6.1 as a possible explanation for the ncf) enhancement. The 
branching ratios are in fair agreement with those from Conforto (1967) and those for 
■k^KlKl, except for the much larger Kir S-wave in tt^K^Kl. 



9 Annihilation at Rest into 5tt 

Resonances in p + p~ have been reported in pn annihilation in deuterium: A 2 ++ state was 
observed in bubble chamber exposures in liquid deuterium (Bridges, 1987). An enhance- 
ment was also seen around 1500 MeV by the Asterix collaboration at LEAR in gaseous 
deuterium but no spin-parity analysis was performed (Weidenauer, 1993). This state 
was interpreted as /2(1565) in its pp decay mode, the slightly lower mass being due to 
7r-rescattering with the recoil proton exciting the A resonance (Kolybashov, 1989). The 
2 ++ assignment was, however, disputed in favour of ++ by a reanalysis of the bubble 
chamber data (Gaspero, 1993). 

The Crystal Barrel has also searched for scalar mesons decaying to Air. The 47r° decay 
mode was investigated using pp annihilation into 5tt°, leading to 10 detected photons 
(Abele, 1996c). The branching ratio for annihilation into 57r° was found to be (7.1zbl.4) x 
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10 -3 . After removing the r\ — ► 3ir° contribution they performed a maximum likelihood 
analysis of a sample of 25,000 5tt° events. The data demand contributions from 7r(1330) — > 
37T° and from two scalar resonances decaying to Att°. The mass and width of the lower 
state, presumably /o(1370), could not be determined precisely. 

The upper state has mass ~ 1505 MeV and width ~ 169 MeV and decays into two 
■kit S-wave pairs and 7r°(1300)7r° with approximately equal rates. The branching ratio 
for pp — > /o(1500)7T°— > 57T° is (9.0 ± 1.4) xlO -4 . Using the 2ir branching ratio from the 
coupled channel analysis one finds 

^0(1500)^4.) _ 2 _ 1±0 _ 6; 



5(70(1500) 2tt) 
5(70(1500) -> 4^°) 
5(70(1500) -» W ) 



1.5 ±0.5. (113) 



The 47T° (2-7T ) contribution in ( |112| ) has been multiplied by 9 (3) to take into account the 
unobserved charged pions. The ratio ( |1 12|) is, in principle, a lower limit which does not 
include pp. However, a reanalysis of the Mark III data on J/ip — ► 727r + 27r~ finds evidence 
for 7o(1500) decaying to 4tt through two S-wave dipions with negligible pp contribution 



(Bugg, 1995). The ratio (113) is in agreement with the result for the former 7o(1590) 
from Aide (1987) in Tr~p -> 4vr n: 0.8 ± 0.3. 

A scalar resonance with mass 1374 ± 38 and width 375 ± 61 MeV decaying to 7r + 7r _ 27r° 
was also reported by Crystal Barrel in the annihilation channel 7r + 7r~37r° (Amsler, 1994d). 
The branching ratio for 7r + 7r~37r° was measured to be (9.7 ± 0.6) %. The Air decay mode 
of the resonance is five times larger than the 27r, indicating a large inelasticity in the 2tt 
channel. The relative decay ratio to pp and two tttt S-waves is 1.6 ± 0.2. However, the 
data do not exclude the admixture of a 7o(1500) contribution. 



10 The New Mesons 

It is instructive to first check the consistency within Crystal Barrel data and also com- 
pare with previous measurements. The squares of the isospin Clebsch-Gordan coefficients 
determine the total branching ratios, including all charge modes. Note that two neutral 
isovectors (e.g. a2(1320)-7r°) do not contribute to I = 1. Table [l^ gives the corresponding 
weights to KK. For K* (892), (1430), a 2 (1320) and a (1450) contributions to KKtt the 
tt°K l K l ratios in Table |l| must be multiplied by 12 while those for 7 2 (1270), 7 (1370), 
7o(1500) and 7 2 (1525) must be multiplied by 4 (since K$ is not observed). For tt^K^K^ 
the resonance contributions in Table [Tl] have to be multiplied by 3 for K*(892), Kq(1430), 
by 3 for a (980), a 2 (1320), /)(1450/1700) from I = and by 2 for a 2 (1320), p(1450/1700) 
from 1=1. 

Furthermore, the branching ratios must be corrected for all decay modes of the in- 
termediate resonances to obtain the two-body branching ratios in Table [l3|. We have 
used the following decay branching ratios (Barnett, 1996): (28.2±0.6)% for 72(1270) — > 
7r°vr , (4.6±0.5)% for 7 2 (1270) KK, (14.5±1.2)% for a 2 (1320) -> r/vr, (4.9±0.8)% for 
a 2 (1320) -> KK, (0.57±0.11)% for a 2 (1320) r/'vr and (88.8±3.1)% for 7 2 (1525) KK. 
There are indications that the a 2 (1320) contribution to ttttt] is somewhat too large. Oth- 
erwise the consistency is in general quite good and Crystal Barrel results also agree with 
previous data. This gives confidence in the following discussion on the new mesons. 
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10.1 a (1450) 

We begin with the isovector ao(1450) which has been observed in its 7771", rj'ir and KK 
decay modes. Averaging mass and width from the coupled channel and the KK analyses 
one finds: 

a (1450) : m = 1474 ± 19, r = 265 ± 13 MeV. (114) 

The ao(1450) decay rates are related by SU(3)-flavor which can be tested with Crystal 
Barrel data. Following Amsler and Close (1996) we shall write for a quarkonium state 

\qq) = cosa|ran) — sina|ss), (115) 

where 



\nn) = (uu 



+ dd)/V2. (116) 



The mixing angle a is related to the usual nonet mixing angle 9 (Barnett, 1996) by the 
relation 

a = 54.7° +0. (117) 

Ideal mixing occurs for 6 = 35.3° (-54.7°) for which the quarkonium state becomes pure 
ss (nn). 

The flavor content of 77 and 7/ are then given by the superposition (see also Eq. (|2l|)) 



\rj) = cos0|nn) — sin0|ss), 

|t? ) = sin0|nn) + coscj)\ss), (118) 

with (p = 54.7° + Op, where 9 P is the pseudoscalar mixing angle which we take as 6 P = 
(-17.3 ± 1.8)° (Amsler, 1992b). 

The partial decay width of a scalar (or tensor) quarkonium into a pair of pseudoscalars 
Mi and M 2 is given by 

T(M 1 ,M 2 )= 1 2 (M l ,M 2 )f L (p)p. (119) 

The couplings 7 can be derived from SU(3)-flavor. The two-body decay momentum is 
denoted by p and the relative angular momentum by L. The form factor 

2 

f L (p)=p 2L eM-{^) (120) 

provides a good fit to the decay branching ratios of the well known ground state 2 ++ 
mesons if (3 is chosen > 0.5 GeV/c (Amsler and Close, 1996). We choose f3 = 0.5 GeV/c 
although the exponential factor can be ignored {(3 — * 00) without altering the forthcoming 
conclusions. Replacing Jl(p)p by prescription ( p8| ) also leads to a good description of decay 
branching ratios provided that pn > 500 MeV/c, corresponding to an interaction radius 
of less than 0.4 fm (Abele, 1997g). 

The decay of quarkonium into a pair of mesons involves the creation of qq pair from 
the vacuum. We shall assume for the ratio of the matrix elements for the creation of ss 
versus uu (or dd) that 

H {0\V\uu) K ' 

This assumption is reasonable since from the measured decay branching ratios of tensor 
mesons one finds p = 0.96 ± 0.04 (Amsler and Close, 1996). Similar conclusions are 
reached by Peters and Klempt (1995). 
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Let us now compare the Crystal Barrel branching ratios for (1320) decays to KK, 
t]tt and tj'tt with predictions from SU(3). One predicts for an isovector with the coupling 
constants 7 given in the appendix of Amsler and Close (1996): 



r(a^(1320) -» K ± K°) 



1 



hipk)p. 



K 



r(a°(1320) - 
r(o^(1320) 



r]TT°) 
^ tj'tt 



2cos 2 <f) h(p v )p v 

f2(PK)PK 



r(a£(1320) -» rj-K 



tan 



f2(p V )Pr) 



0.295 ±0.013, 



0.029 ± 0.004. 



(122) 
(123) 



Both ratios are in excellent agreement with world data (Barnett, 1996). The Crystal 
Barrel numbers to be compared with are taken from Tables [H] and 11 for 1 = 0. One 
finds the ratios of branching ratios 



B(af (1320) -> K±K°) 
5(a§(1320) -»■ r]7T°) 



5(a§(1320) 



rj 7r 



B(a§(1320) 



21 +a04 
u - zi -0.06> 



0.034 ± 0.009, 



r/7r u 



(124) 
(125) 



in fair agreement with SU(3) predictions. We now compare predictions and data for 
ao(1450). From SU(3) one expects, using Eq. ( ggg ) and ( |T23D with L = 0: 

T(a^(1450) -» K±K°) 



r(ag(1450) - 
r(ag(1450) 



r]ir v ) 
-> tj'tt ) 



T(ag(1450) -> ??^ ) 



0.72 ±0.03, 
0.43 ± 0.06, 



in agreement with the experimental results from Tables |10| and |11] : 

5(^(1450) -» ^i^T ) 



B(ag(1450) r/vr c 
S(ag(1450) 



r/'vr ) 



5(ag(1450) 7?vr c 



0.87 ±0.23, 



0.34 ±0.15. 



(126) 
(127) 

(128) 
(129) 



This then establishes ao(1450) as a qq isovector in the scalar nonet. 

The existence of ao(1450) adds further evidence for ao(980) being a non-qq state. The 
/o(980) and ao(980) have been assumed to be KK molecules (Weinstein, 1990; Close, 
1993). This is motivated by their strong couplings to KK - in spite of their masses close 
to the KK threshold - and their small 77 partial widths. For /q(980) the 27 partial width 
is T 77 = (0.56 ± 0.11) keV (Barnett, 1996). The relative ratio for a (980) decay to KK 
and r]7t has been determined by Crystal Barrel (Eq. (|111|)). Using 



r 77 5(a (980) 



r]Tr) = 0.24 ± 0.08 keV 



(130) 



(Barnett, 1996), one then derives the partial width T 77 = (0.30 ± 0.10) keV. Thus the 27 
partial widths for /o(980) and ao(980) appear to be similar, close to predictions for KK 
molecules (0.6 keV) and much smaller than for qq states (Barnes, 1985). However, not 
everybody agrees: In Tornqvist (1995) the / (980)// (1370) and the a (980)/a (1450) 
are different manifestations of the same uniterized ss and ud states, while the broad 
structure around 1100 MeV is the uu + dd state (Tornqvist and Roos, 1996). This then 
leaves _/b(1500) as an extra state. 
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10.2 / (1370) and / (1500) 

From the single channel analyses and the KK decay mode we find for /q(1370) and 
/o(1500) the average masses and widths: 

/o(1370) : m = 1360 ± 23 MeV, T = 351 ± 41 MeV, 

/o(1500) : m = 1505 ± 9 MeV, T = 111 ± 12 MeV. (131) 

The closeness of ao(1450) and /o(1500) or even / (1370) masses is conspicuous and points 
to a close to ideally mixed scalar nonet, one of the latter mesons being one of the qq 
isoscalars. However, /o(1500) with a width of about 100 MeV is much narrower than 
a (1450), / (1370) and ^(1430) with widths of typically 300 MeV. Theoretical predic- 
tions for the widths of scalar qq mesons, based on the 3 Po model, agree that scalar qq 
mesons have widths of at least 250 MeV (for a discussion and references see Amsler and 
Close (1996)). We therefore tentatively assign /q(1370) to the ground state scalar nonet. 

If /o(980) is indeed a molecule then the (mainly) ss member of the scalar nonet still 
needs to be identified. We now show from their decay branching ratios that neither 
/o(1370) nor /o(1500) are likely candidates. To investigate the quark content of /o(1500) 
we calculate its relative couplings to 7777, 7777' and KK and search for a common value of 
the scalar mixing angle a. The ratios of couplings for a pseudoscalar mixing angle cj) are 
(Amsler and Close, 1996): 

Ri 

R2 



7 2 (w) (cos 2 (j) — v2 tana sin 2 </> 

7 2 (7T7r) 3 

7 2 (?7?/) 2(cos0sin0[l + \/2tanaP 2 

7 2 (7T7r) 3 



j 2 (KK) (1 - v^tan 



o 



,2 



7T7T : 


(29.0 


± 


7.5) 
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7777 : 


(4.6 
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1.3) 


% 


7777' : 


(1.2 


± 


0.3) 
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KK : 


(3.5 


± 


0.3) 


% 


47r : 


(61.7 


± 


9.6) 


% 



^3 = ' h r = -• (132) 

For 7777 and ttit we use the branching ratios from the coupled channel analysis (Table |10| ) 
and multiply irn by 3 to take into account the 7r + 7r _ decay mode. The branching ratio for 
KK is taken from Table O and is multiplied by 4. The branching fractions are, including 
the 47r mode from Eq. (|112j) and ignoring a possible small pp contribution to Air: 



(133) 



After correcting for phase space and form factor (Eq. ( |120[ )) we obtain: 

Ri = 0.195 ± 0.075, R 2 = 0.320 ± 0.114, R 3 = 0.138 ± 0.038. (134) 

Since /o(1500) lies at the 7777' threshold we have divided the branching ratios by the phase 
space factor p integrated over the resonance and have neglected the form factor when 
calculating i?2. 

Previously, the upper limit for R3 was < 0.1 (95% confidence level) from Gray (1983), 
in which case no value for the mixing angle a could simultaneously fit R\, R2 and R3 
(Amsler and Close, 1996), therefore excluding / (1500) as a qq state. The (la) allowed 
regions of tana are shown in Fig. ^ for the ratios (|134| ). The agreement between R± and 
R3 is not particularly good. Remember, however, that branching ratios are sensitive to 
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interference effects and therefore caution should be exercized in not overinterpreting the 



apparent discrepancy in Fig. 34. On the basis of the ratios ( |134| ), one may conclude that 
/o(1500) is not incompatible with a mainly uu + dd meson (a = 0). For a pure ss state 
(a = 90°) the ratios ( |134| ) would, however, become infinite. Therefore /o(1500) is not the 
missing ss scalar meson. 

Similar conclusions can reached for /o(1370) which has small decay branching ratios to 
nn and KK. Precise ratios Ri are difficult to obtain in this case since the branching ratio 
to tttt in Table |l0| also includes the low energy tttt S-wave, in particular /o(400 — 1200). 

This analysis shows that both /o(1370) and /o(1500) are compatible with isoscalar 
uu + dd states, although the latter is much too narrow for the ground state scalar nonet. 
This then raises the question on whether /o(1500) could not be the first radial excitation 
of /o(1370). This is unlikely because (i) the splitting between ground state and first radial 
is expected to be around 700 MeV (Godfrey and Isgur, 1985), (ii) the next Kq lies at 1950 
MeV (Barnett, 1996) and, last but not least, first radials are expected to be quite broad 
(Barnes, 1997). 

The most natural explanation is that /o(1500) is the ground state glueball predicted 
in this mass range by lattice gauge theories. However, a pure glueball should decay to 
7T7T, 7/r/, 7/7/' and KK with relative ratios 3 : 1 : : 4, in contradiction with our ratios 
Ri. In the model of Amsler and Close (1996) the finite nn' and the small KK rates can 
be accommodated by mixing the pure glueball Go with the nearby two nn and sH states. 
Conversely, the two isoscalars in the qq nonet acquire a gluonic admixture. In first order 
perturbation one findsf ! 

\ffiKM\W \Go) + Z(V2\nn) + u\ss}) 

l/o(1500)) = TTTHfT^) ' (135) 

where oj is the ratio of mass splittings 

m(G )-m(nn) 

u= !n x T—r- (136) 

m(Go) — m{ss) 

In the flux tube simulation of lattice QCD the pure gluonium Go does not decay to irir nor 
to KK in first order and hence /o(1500) decays to tttt and KK through its qq admixture 
in the wave function. If Go lies between the two qq states, u is negative and the decay to 
KK is hindered by negative interference between the decay amplitudes of the nn and sis 
components in Eq. (|135| ), The ratio of couplings to KK and tttt is 



-f\KK) = (1 + 

7 2 (7T7r) 3 



00 



,2 



(137) 



The cancellation is perfect whenever Go lies exactly between nn and ss (w = —1). We 
find with the measured R3 two solutions, uj = - 0.36 or -1.64. Assuming that /o(1370) 
is essentially nn (with a small gluonic admixture) this leads to an ss state around 1900 
MeV or 1600 MeV, respectively. Furthermore, the ratio of tttt partial widths for /o(1500) 
and /o(1370), divided by phase-space and form factor, is given by 

f 3,1/0(1500)1 = 0Sj pa) 

r 2 „[/o(1370)l l+fp + u*) 



using for /q(1500) the branching ratios from Tables [LC] and 11 and the 4tt/2tt ratio (|TT 



There is, however, a large uncertainty in the ratio ( 138 ) due to the branching ratios of 



10 We assume here that the quark-gluon coupling is flavor blind, see Amsler and Close (1996) for a 
generalization. 
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/o(1370) which cannot easily be disentangled from /o(400 — 1200). This then leads to 
|£| ~ 0.6 and according to Eq. (135) to about 30% or about 60% glue in /o(1500) for an 
ss state at 1600 MeV or 1900 MeV, respectively. 

The /o(1500) has also been observed in pp annihilation at higher energies (Armstrong, 
1993) and in other reactions, in particular in central production, decaying to 2ir + 2ir~ 
(Antinori, 1995). The VES experiment, studying ir~p interactions on nuclei at 36 GeV/c, 
has reported a resonance, 7r(1800), decaying to ir~r]r]' (Beladidze, 1992) and ir~r]r] (Amelin, 
1996). The 7r(1800) appears to decay into a resonance with mass 1460 ± 20 MeV and 
width 100 ± 30 MeV - in agreement with /o(1500) - with a recoiling ir. They report an 
rjrj'/rjri ratio of 0.29 ± 0.07 which is in excellent agreement with the Crystal Barrel ratio 
for / (1500) decays, 0.27 ± 0.10 (Table 0). Note that if vr(1800) is indeed a qqg (hybrid), 
as advocated by Close and Page (1995), then decay into gluonium is favoured,. 

A reanalysis of J/ip radiative decay to 2tt + 2tt~ finds evidence for /o(1500) decaying 
to two S-wave dipions with a branching ratio in J/ip — > jAtt of (5.7 ± 0.8) x 10~ 4 (Bugg, 
1995). This leads to an expected branching ratio of (2.7±0.9) x 10 -4 in J/ip — > ^2ir, using 
the Crystal Barrel result ( |112| ). It is interesting to compare this prediction with data on 
J/ip —* 77r + 7T~ from Mark III (Baltrusaitis, 1987) where /2(1270) is observed together 
with a small accumulation of events in the 1500 MeV region. Assuming that these are 
due to /o(1500), one finds by scaling to /2(1270) a branching ratio in J/ip — > ^/2tt of 
~ 2.9 x 10~ 4 , in agreement with the above prediction. 

Summarizing, /o(1500) has been observed in pp annihilation in several decay modes, 
some with very high statististics (~ 150,000 decays into ) and also in other processes 
that are traditionally believed to enhance gluonium production, central production and 
J/ip radiative decay. The -Kq(1430) and oo(1450) define the mass scale of the qq scalar 
nonet. The /o(1500) is not the missing sH and is anyway too narrow for a scalar qq state. 
The most natural explanation for /o(1500) is the ground state glueball mixed with nearby 
scalars. The missing element in this jigsaw puzzle is the ss scalar expected between 1600 
and 2000 MeV. The analysis of in flight annihilation data will hopefully provide more 
information in this mass range. The spin of /j(1710) has not been determined unambigu- 
ously. If J = is confirmed then /j(1710) could be this state or, alternatively, become 
a challenger for the ground state glueball (Sexton, 1995). A more detailed discussion on 
/o(1500) and /j(1710) can be found in Close (1997). 



10.3 / 2 (1565) 

The /2(1565) with mass 1565 ± 20 and width 170 ± 40 MeV has been observed first by 
the Asterix collaboration at LEAR in the final state 7r + 7r~7r° in hydrogen gas (May, 1989, 
1990) and then by Aker (1991) in the 3ir° final state in liquid. The 3ir° analysis gave 60% 
P-state contribution to 3tt° in liquid with roughly equal intensities from /2(1565) of 9%, 
each from 1 Sq, 3 Pi and 3 P2- The full 37r° data sample demands 46% P-state annihilation 
(Amsler, 1995f) and also requires a tensor around 1530 MeV (section p.3[ ). The ss tensor, 
^(1525), has been observed in its KK decay mode (section |8.7|) . From the observed rate 
(Table |l3| ) and the known ^(1525) decay branching ratios (Barnett, 1996) one finds that 
,/2(1525) cannot account for much of the 2 ++ signal in tttt or 7777. Hence /2(1565) is not 
$(1525). 

The large P-state fraction in the 3ir° channel in liquid is not too surprising: the 
corresponding channel pp — > 7r + 7r~7r in liquid proceeds mainly from the (I =0) 3 Si atomic 
state while the (7=1) 1 So is suppressed by an order of magnitude (Foster, 1968b), as are 
normally P-waves in liquid. On the other hand, the channel pp — > 3ir° proceeds only 
through the (7=1) 1 Sq atomic state while 3 Si is forbidden. Indeed the branching ratios 
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for 37r° is an order of magnitude smaller than for ir + ir~ir° (Table ^). Hence for 3ir° S- 
and P-wave annihilations compete in liquid. 

A fraction of 50% P-wave was also required in the Dalitz plot analysis of the 1=1 
final state 7r~7r°7r° at rest in liquid deuterium (Abele, 1997h) which shows evidence for 
/o(1500) and /2(1565) production and requires in addition the p-meson and two of its 
excitations, p _ (1450) and p~(1700), decaying to ir^ir . 

However, the coupled channel analysis described in section |8.4| , ignoring P-waves, still 
requires a tensor at 1552 MeV. Neglecting P-waves increases slightly the contribution 
from /o(1500) while decreasing the contribution from /2(1565), although the rates remain 
within errors (compare the two 7r°7r° branching ratios in Table |l^ for the single and 
coupled channel analyses). 

The alternative N/D analysis also reproduces the features of the 3ir° Dalitz plot with- 
out P-wave contributions, in particular the scalar state around 1500 MeV (Anisovich, 

1994) . A tensor contribution with mass ~ 1565 and width ~ 165 MeV is, however, still 
required, but most of the blob structure in Fig. |2l](c) is taken into account by interferences 
in the low energy irir S-wave. 

An N/D analysis of Crystal Barrel data, together with former data from other reac- 
tions, also reports a tensor with mass 1534 db 20 and width 180 db 60 MeV (Abele, 1996a). 
They report strong pp and ww contributions and therefore assign this signal to /2(1640) 
discovered by GAMS in ir~p — ► u>un (Aide 1990), also reported to decay into Air by the 
Obelix collaboration in np — > 5ir (Adamo, 1992). It should be emphasized, however, that 
in Abele (1996a) the inelasticity in the K-matrix is attributed to pp and ww, although no 
47r data are actually included in the fit. Given that mass and width of the tensor agree 
with /2(1565), but disagree with /2(1640), it seems more natural to assign the 2 ++ signal 
to the former. It is interesting to note that a 2 ++ pp molecule (Tornqvist, 1991) or a 2 ++ 
baryonium state (Dover, 1986) decaying strongly to pp (Dover, 1991) are predicted in this 
mass range. The /2(1565) could be one of these states. 

In conclusion, there is a certain amount of model dependence when extracting the 
precise production and decay rates of /2(1565). However, annihilation data require both 
a scalar and a tensor around 1500 MeV and the parameters and rates for /o(1500) are 
reasonably stable, independent of /2(1565) contribution. The Crystal Barrel data in gas 
will hopefully settle the issue of the fraction of P-wave in three-body annihilation. 

11 E/i Decay to r]nn 

The E meson, a ^ state, was discovered in the sixties in the KKir mass spectrum of 
pp annihilation at rest into (KsK^ir^)ir + ir~ . Its mass and width were determined to 
be 1425 d= 7 and 80 ± 10 MeV (Baillon, 1967). Its quantum numbers have remained 
controversial since other groups have claimed a ^ state (now called 7/(1440)) and a 1 ++ 
state (now called /x(1420)) in this mass region from various hadronic reactions. A broad 
structure (previously called t), has also been observed in radiative J/ijj decay to KKir 
(Scharre, 1980). Initially determined to be '", the E/i structure was then found to 
split into three states, the first (0 '") at 1416 ± 10 MeV decaying to ao(980)7r, the second 
(presumably the 1++ /i(1420)) at 1443 ± 8 MeV and the third (0 _+ ) at 1490 db 16 MeV, 
both decaying to K*K (Bai, 1990). The widths were not determined accurately. The 
Obelix collaboration has analyzed the KKir mass spectrum in pp annihilation at rest in 
liquid and has also reported a splitting of the E meson into two pseudoscalar states at 
1416 db 2 MeV (r = 50 ± 4 MeV) and 1460 ± 10 MeV (r = 105 db 15 MeV) (Bertin, 

1995) . We shall refer to these pseudoscalars as 7/(1410) and r/(1460). When using gaseous 
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hydrogen, one expects the production of 1 ++ mesons from 3 Pi recoiling against an S- 
wave dipion: Bertin (1997) indeed observes three states in the KKtt mass spectrum in 
gas: 7/(1410), /i(1420) and 77(1460). 

In J/ip radiative decay, t decays to KKtt through the intermediate ao(980) and hence 
a signal was also expected in the ao(980)7r — * t]tttt mass spectrum. This has indeed been 
observed by Mark III and DM2: Bolton (1992) reports a signal in clq tt^ at 1400 ± 6 MeV 
(r = 47 ± 13 MeV) and Augustin (1990) in r/vr+vr- at 1398 ± 6 MeV (T = 53 ± 11 MeV). 
We shall tentatively assign these signals to 7/(1410). 

To clarify whether the structures observed in J/i/j radiative decay and in pp annihi- 
lation are compatible and in particular to confirm the quantum numbers of E (0 ^ and 
not 1 ++ ) Crystal Barrel has searched for the t]tttt decay mode of the E meson in the 
reaction pp — > (7]tt + tt^)tt tt and (?77r 7r )7r + 7r~, leading to two charged particles and 67 
(Amsler, 1995h; Urner, 1995). Since the rate for this reaction was expected to be rather 
small (~ 10~ 3 of all annihilations), an online trigger required 8 clusters in the barrel and 
at least two ir° and one 77. A 7C kinematic fit then selected the channel 7r + 7r~27r r/ while 
suppressing tt + Tt~ 2r]TT° and 7r + 7r~37r°. The branching ratio for 7r + 7r~27r°77 was found to 
be (2.09 ± 0.36) %. 

The final state tt + tt~2tt°7] includes a strong contribution from ujr]7r (u) — ► 7r + 7r~7r°), a 
channel that has been studied in its 77 decay mode (Amsler, 1994c). Events compatible 
with ujtjtt were removed, leaving a sample of about 127,000 7r + 7r~27r°77 events. The 
evidence for 77(1410) decaying to t]tttt is shown in the TT°TT°r] and ir + TT~rj mass distributions 
(Fig. ^). Some 9,000 7/(1410) decays into r/TTir are observed in Fig. |3l| an order of 
magnitude more than for E to KKtt in the seminal work of Baillon (1967). 

A partial wave analysis was performed using a maximum likelihood optimization. The 
signal at 1400 MeV was described by the annihilation channels 



pp -> f7(1410)(-> 
-> 77(1410)(-> 
-> 77(1410)(— s- 
-> fj(1410)(-> 

where cr° and are shorthands for the 7r°7r° and tt + tt~ S-waves. The latter were 



described by prescriptions of the form (100) which is reasonable since the txtt masses lie 
below 900 MeV. Background contributions, e.g. from rjp°a° and r]'p°, were also included 
in the fit. Figure 36 shows for example the a^(980) angular distribution in the r/(1410) 
rest frame together with the best fit for a ^ state. The data exclude 1 ++ , hence r/(1410) 
is definitively pseudoscalar and is produced from the 1 Sq atomic state. It has mass and 
width 

7/(1410) : m = 1409 ± 3 MeV, r = 86 ± 10 MeV. (140) 

The width is somewhat larger than for tjtttt in J ftp decay (Bolton, 1992; Augustin, 1990) 
and for KKtt in pp annihilation at rest (Bertin, 1995). 

The branching ratio to KKtt has been measured earlier (Baillon, 1967): 

B(pp -» Etttt, E -» KKtt) = (2.0 ± 0.2) x 10" 3 , (141) 



while Crystal Barrel finds for the t]tttt mode 



ag(980)vr )(j + -, 

0^(980)^)0-°, (139) 



B(pp -> 77(1410)vr7r, 7/(1410) t/vttt) = (3.3 ± 1.0) x 10" 3 . (142) 

The fit yields rja and ao(980)(^ t]tt)tt decay contributions with a relative rate of 0.78 
± 0.16. Assuming that 50% of the KKtt mode proceeds through 7/(1410) decaying to 
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ao(980)-7r (Baillon, 1967) one can estimate from these branching ratios 



B(a 980 -> KK) , , 

— ^ ^~ 0.54 ±0.18, 143) 

5(a (980) r?vr) ' v ; 

somewhat larger but not in violent disagreement with the result ( |111| ). 

The observation of the r/7r°7r° decay mode also lifts the earlier isospin ambiguity for 
the E meson and clearly establishes that this state is isoscalar (C = +1). Note that the 
Crystal Barrel data do not exclude the presence of the other 1 = pseudoscalar, 77(1460), 
since the latter was observed in K*K and not in rjTrir. 

The first radial excitation of the r] could be 77(1295) decaying to rjirn (Barnett, 1996). 
Hence one of two pseudoscalars in the i structure could be the radial excitation of the 
r/. The near equality of the 7/(1295) and 7r(1300) masses suggests an ideally mixed nonet 
of h radials. This implies that the second isoscalar in the nonet should be mainly ss 
and hence decays to K*K, in accord with observations for 77(1460). This scheme then 
favours an exotic interpretation for 77(1410), perhaps gluonium mixed with qq (Close, 1997) 
or a bound state of gluinos (Farrar, 1996). The gluonium interpretation is, however, 
not favoured by lattice gauge theories, which predict the state above 2 GeV (see 
Szczpaniak (1996)). 



12 Summary and Outlook 

Crystal Barrel has collected 10 8 pp annihilation at rest in liquid hydrogen, three orders of 
magnitudes more than previous bubble chamber experiments. The results reviewed in this 
report were achieved thanks to the availability of pure, cooled and intensive low energy 
antiproton beams which allow a good spatial definition of the annihilation source and 
thanks to the refinement of the analysis tools warranted by the huge statistical samples. 
The data processed so far concentrate on annihilations at rest into 0-prong that had 
not been investigated before. The data collected with the 0-prong trigger correspond to 
6.3 x 10 8 annihilations. 

The measurement of the branching ratio for annihilation into 7T 7T° leads, together 
with a cascade calculation of the antiprotonic atom, to a fraction of (13 ± 4) % P- wave 
in liquid hydrogen. Therefore S-wave dominance has been, in general, assumed to analyze 
the data. 

The branching ratios for annihilation into two neutral light mesons (vr°77, 7r°7/, 7777, 7777', 
lot], air)', r)p°, r]'p°) reveal the interplay of constituent quarks in hadrons. The non-planar 
quark rearrangement graph must play an important role in the annihilation process. Using 
the OZI rule the pseudoscalar mixing angle was determined to be (—17.3 ± 1.8)°. 

However, the production of 4> mesons is enhanced in nearly all channels compared to 
predictions by the OZI rule. The most significant deviation is found in the annihilation 
channel ir°(p. After phase space correction, the k^^/it^uj ratio is (10.6 ± 1.2) % while OZI 
predicts 0.42 %. Whether this enhancement can be explained by final state corrections or 
by ss pairs in the nucleon is not clear yet. The analysis of data in P-state annihilations, 
in gas or at higher momenta, will be helpful in settling the nature of this phenomenon. 

In electromagnetic processes, the radiative annihilations vr°7, 777 and toj have been 
observed with rates consistent with predictions from VDM, but (fry may be enhanced. 
The branching ratio for u — > 777, (6.6 db 1.7) x 10 -4 , was measured independently of p — u 
interference. This results solves the ambiguity in e + e~ formation experiments, selecting 
the constructive p — u interference solution. The 77 — » 37r° Dalitz plot is not homogeneous 
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but shows a negative slope of a = 0.052 ± 0.020. Crystal Barrel data also confirms the 
evidence for the direct decay r/ — > 7r + 7r~7, in addition to if — ► /ry. 

A decisive progress has been achieved in understanding scalar mesons by studying 
annihilation into three pseudoscalars. An isovector state, ao(1450), with mass and width 
(m, T) = (1474 ± 19, 265 ± 13) MeV has been observed to decay into i]7r, tj'tv and KK 
with rates compatible with SU(3) flavor. The existence of ao(1450) adds evidence for 
ao(980) not being qq, but perhaps a KK molecule. The ratio of ao(980) decay rates to 
KK and rjir was measured to be 0.24 ± 0.06. 

An isoscalar state, /o(1370), with mass and width (m, T) = (1360 ± 23, 351 ± 41) MeV 
has been observed to decay into tttt, r/r), KK and Air. Obtaining accurate branching ratios 
for Jo (1370) is difficult due to interferences with the broad structure /o(400 — 1200). The 
nature of this structure (meson or slowly moving background phase) is unclear. Whether 
/o(400 — 1200) is really distinct from /o(1370) is not entirely clear. Both questions will 
probably remain with us for some time. The states ao(1450), /o(1370) and Xq(1430) are 
broad, consistent with expectations for qq scalar mesons. The small coupling of /o(1370) 
to KK makes it an unlikely candidate for the ss meson, which is therefore still missing. 
More information on this state will hopefully emerge in the mass range above 1600 MeV 
from Crystal Barrel data in flight. 

An additional isoscalar state, /o(1500), with mass and width (m, V) = (1505 ± 9, 
111 ± 12) MeV has been observed to decay into tttt, rjrj, rjrf, KK and 4tt. The decay 
branching ratios are 29, 5, 1 , 3 and 62 %, respectively. These rates exclude this state to 
be the missing ss. Hence /o(1500) is supernumerary and anyway too narrow to be easily 
accommodated in the scalar nonet. The likely explanation is that /o(1500) is the ground 
state glueball predicted by QCD, mixed with the two nearby qq isoscalars, /o(1370) and 
the higher lying sH state. 

The tensor /2(1565) is dominantly produced from P-states. It is, however, still required 
to fit the data when assuming pure S-wave annihilation: Its mass and width are ~ 1552 
and ~ 142 MeV. The systematic inclusion of P-wave annihilation at rest in all analyses 
is, however, prevented by the large number of fit parameters. Data from Crystal Barrel 
and Obelix in liquid and gaseous hydrogen might alleviate this problem, perhaps also 
modifying slightly some of the branching ratios obtained from liquid only. 

A 0~ + state, r?(1410), with mass and width (m,T) = (1409 ± 3, 80 ± 10) MeV has 
been observed to decay into r]-K°iT and t]tt + tt~ with approximately equal rates through 
ao(980)-7r and (irir)sr] . The neutral decay mode establishes this state as an isoscalar and, 
together with other experiments, strengthens the evidence for two 7 = pseudoscalars in 
the 1400 - 1500 MeV region. 

Crystal Barrel data also show evidence for the two excitations of the p meson, /?(1450) 
and /?(1700) and for the radial excitation 0^(1650). The analysis of in-flight data will 
hopefully reveal further radial excitations, hybrid mesons and higher mass glueballs. For 
glueballs, a more definitive progress will probably be achieved in radiative J/tp decay 
at a high luminosity e + e~ factory or in central collisions at the forthcoming Compass 
experiment at CERN. 
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Table 1: Summary of data (in millions of events) with a minimum bias trigger (MB), for 0-, 2-, 
4-prong and with more specialized triggers at rest (first three rows) and in liquid hydrogen at 
high p momenta (in MeV/c). 

LH 2 : liquid hydrogen; LD 2 : liquid deuterium; GH 2 : gaseous hydrogen. 
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Table 2: Branching ratios B for pp annihilation at rest in liquid. See Amsler and Myhrer (1991) 
for annihilation in gaseous hydrogen. Further branching ratios from Dalitz plot analyses are 



listed in Table [13] below. 
a From ui — > 7r°7 
6 From oj — > 7T + 7r _ 7r° 

c average between Baltay (1966), Espigat (1972) and Foster (1968a) 
\ Crystal Barrel experiment 
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Table 3: Pseudoscalar mixing angle 6 P derived from the measured ratios of two-body branching 
ratios (0j = 35.3°). The first four rows assume only the QLR in the annihilation process. The 
last six rows assume in addition the dominance of the annihilation graph A. 

Ratio Prediction 9„ \°] 
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Table 4: Damping factors F L {p) where z stands for (p/pn) 2 and p R is usually taken as 197 
MeV/c (after Hippel and Quigg (1972)). 

L F L (p) 







I / 22 



2 + 1 



2 y ^ 



(2-3) 2 +92 



o / 2772~ 

° V 2(2-15) 2 +9(22-5) 2 

4 / .. 12,746z 4 

V (2 2 -452+105) 2 +252(22-21) 2 



Table 5: Branching ratios -B for radiative pp annihilation at rest in liquid from Crystal Barrel 
(Amsler, 1993c, 1995c). The lower and upper limits L and U, calculated from VDM, are given 
in the third and fourth column, respectively, 
f 95% confidence upper limit. 
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Table 6: Branching ratios for production at rest in liquid. 
a updates Amsler (1995c) 

b annihilation in gas extrapolated to pure S-wave annihilation 
d using Chiba (1988) in liquid 
d using Bizzarri (1971) in liquid 
\ Crystal Barrel experiment 
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Table 7: Ratio of to u production in low energy annihilation in liquid. 
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Table 8: Relative sign of a± and a 2 for pp annihilation into nKK (see text). 

Channel 1 S G {pp) 3 <Si(pp) 

1=0 1=1 I=Q 1=1 

^K^K + 1 1 + 

7T K+K- + + 
Tr°K°KV + + 
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Table 9: Branching ratios for pp annihilation at rest into three narrow mesons. Mesons in 

parentheses were not detected. 

a using B(n°(j)) from Table § and Eq. (g7p 

b average between Armenteros (1965) and Barash (1965) 

f Crystal Barrel experiment 
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Table 10: Branching ratios for pp annihilation at rest in liquid determined from Dalitz plot 
analyses. The branching ratios include all decay modes of the final state stable particles (e.g. 
7T°, rj, 77') but only the decay mode of the intermediate resonance leading to the observed final 
state. 



Channel 


Contributing resonances 


Subchannel 


Branching ratio 


7r u r/r/ 


a (980), / (1370), / (1500), a 2 (1320), X 2 ( 1494) 


/o(1370)vr u 


- 3.5 x 10" 4 


/o(1500)7T° 


(5.5 ± 1.3) x 10~ 4 


02(1320)77 


~ 5.6 x 10~ 5 


X 2 (1494)7r° 


~ 4.0 x 10~ 4 (dominantly P-wave annihilation) 


ttVt? 


a (980), a (1450), a 2 (1320), a' 2 (1650), (j}tt) p 




{tttt) s = /o(400 - 1200) + / (980) + / (1370) 


a (980)7T u 


(8.7 ± 1.6) x 10" 4 


a (1450)7T 


(3.4 ±0.6) x 10~ 4 


(7T7r) S 7] 


(3.4 ±0.6) x 10~ 3 


(lJ7r)p TT 


~ 1.0 x 10~ 4 


a 2 (1320)7r° 


(1.9 ±0.3) x 10~ 3 


a 2 (1650)7r° 


~ 1.3 x 10 -4 


TT TT TT 


(")s + /o(1500), / 2 (1270), / 2 (1565) 


(tttt)s 7T° 


~ 2.6 x 10- 3 


/o(1500)7T° 


(8.1 ±2.8) x 10~ 4 


/ 2 (1270)tt 


~ 1.8 x 10~ 3 


/ 2 (1565)^° 


~ 1.1 x 10" 3 


vr u r/r/, 7r u 7r u r/, 3vr u 


Coupled channels (S-wave annihilation only) 




(7777) s = /o(400 - 1200) + /o(1370) 


(7T7r)5 7T° 


(3.48 ± 0.89) x 10" 3 


(vrTr)^ 7? 


(3.33 ± 0.65) x 10~ 3 


(rir/)^ 7T° 


(1.03 ±0.29) x 10~ 3 


/o(1500)(- tt tt )tt 


(1.27 ±0.33) x 10~ 3 


/o(1500)(-> w)tt° 


(0.60 ±0.17) x 10~ 3 


/ 2 (1270)(-> 7T 7r )7T 


(0.86 ± 0.30) x 10~ 3 


/ 2 (1565)(-> 7r°7r°)7r 


(0.60 ± 0.20) x 10~ 3 


/ 2 (1565)(-> 7777)7^ 


(8.60 ± 3.60) x 10~ 5 (may be more than one object) 


/ r\r~ \r\\ / \ 

a (980)(^ 7r u 77)7r u 


(0.81 ± 0.20) x 10~ 3 


a (980)H 7r°77)r7 


(0.19 ±0.06) x 10~ 3 


a (1450)(^ ^°77)7r° 


(0.29 ±0.11) x 10~ 3 


a 2 (1320)H A)tt° 


(2.05 ±0.40) x 10~ 3 (including a' 2 (1650)) 


TT^rjrj' 


/o(1500) 


/o(1500)vr u 


(1.6 ±0.4) x 10" 4 


7T°7rV 


(7T7r) 5 , a 2 ( 1320), a ( 1450) 


(vrvr) s rf 


(3.1 ±0.4) x 10" 3 


a 2 (1320)^° 


(6.4 ± 1.3) x 10" 5 


a (1450)7r° 


(1.16 ±0.47) x 10~ 4 
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Table 11: Branching ratios for pp annihilation at rest in liquid into kaonic channels. The 
branching ratios include only the decay mode of the intermediate resonance leading to the 
observed final state. 

a from the corresponding K s channels (Armenteros, 1965; Barash, 1965) 
b includes low energy Kn scattering 
c fixed by n ± K*K L data 



Channel 


pp(I) 


Contributing resonances 


Subchannel 




Branching ratio 


^K L K L 
Total" 




K*(892), 7^(1430), a 2 (1320),/ 2 (1270), / 2 (1525) 
fo (1370), f (1500), a (1450) 
(7.5 ± 0.3) x 10" 4 


K*(892)K 

K^{Um)K 

a 2 (1320)7r° 

a (1450)7r° 

/ 2 (1270)tt° 

/^(1525)7r° 

/o(1370)vr 

/o(1500)vr° 


%(0,1) 

%(0) 

^o(l) 


(8.71 ±0.68) x 10~ 5 
(4.59 ± 0.46) x 10~ 5 6 
(6.35 ±0.74) x 10~ 5 
(7.35 ± 1.42) x 10~ 5 c 
(4.25 ± 0.59) x 10~ 5 
(1.67 ±0.26) x 10~ 5 
(2.20 ±0.33) x 10~ 4 
(1.13 ±0.09) x 10~ 4 


Total" 




K* (892), ^(1430), a 2 ( 1320) 
a (980), a (1450), ^(1450/1700) 
(2.73 ±0.10) x 10" 3 


K*(892)K 

K$(1430)K 

a (980)vr 

a 2 (1320)7r 

a (1450)7r 

K*(892)K 

K$(1430)K 

/9(1450/1700)7r 


^o(O) 
%(!) 


(2.05 ±0.28) x 10~ 4 
(8.27 ± 1.93) x 10~ 4 6 
(1.97 1 ° f 4 ) x 10- 4 
(3.99 1 x 10- 4 
(2.95 ± 0.56) x 10- 4 
(3.00 ± 1.10) x 10~ 5 
(1.28 ±0.55) x 10~ 4 b 
(8.73 t lf 5 ) x 10- 5 


K*(892)K 
p(1450/1700)7T 
K*{892)K 
a 2 (1320)7r 


a Si(0) 
3 ^i(l) 


(1.50 ±0.41) x 10" 4 
(8.73 ± 2.75) x 10~ 5 
(5.52 ±0.84) x 10~ 4 
(1.42 ± 0.44) x 10~ 4 



Table 12: Weights of the channels contributing to pp — > irKK. I refers to the pp isospin and i 
to the KK isospin. 

Channel i = 1/2 i = 1 i = i = 1 

7 = 0,1 I=l I=\ I = Q 

^K^K 4 2 4 

it°K + K- 1 11 

-k°K°W 1 11 
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Table 13: Branching ratios B for two-body pp annihilation at rest in liquid hydrogen (including 

all decay modes), calculated from the final states given in the last column. 

a assumes 100% 61 decay to -koj 

b using B(n (j)) from Table § and Eq. (J4Tf) 



Channel 






B 






Final state or ref. 


/ 2 (1270)tt u 


3.1 


± 


1.1 


10" 


<> 

-0 


n t\ t\ 

7T U 7r U 7T U 


/ 2 (1270)vr° 


3.7 


± 


0.7 


10- 


-3 


^K L K L 




4.3 


± 


1.2 


10- 


q 
-0 


Foster (1968b) 


/ 2 (1270)w 


3.26 


± 


0.33 


10- 


n 
- z 


Bizzarri (1969) 


/ 2 (1270)w 


2.01 


± 


0.25 


10- 


_2 


Amsler (1993a) 


/o(1500)7r° 


1.29 


± 


0.11 


10" 


-2 


tt K l K l 


/ 2 (1525)tt 


7.52 


± 


1.20 


10" 


-5 


7T°K L K L 


a 2 (1320)?r 


3.93 


± 


0.70 


10- 


-2 


ttVt? 




3.36 


± 


0.94 


10- 


-2 


ttVt/ 




1.55 


± 


0.31 


10- 


-2 


7T°K L K L 




2.44 


+ 


0.44 

0.64 


10- 


-2 






1.32 


± 


0.37 


10- 


-2 


Conforto (1967) 


0^(1320)^ 


5.79 


± 


2.02 


10" 


-3 




( 3 Si) 


4.49 


± 


1.83 


10- 


-3 


Conforto (1967) 


bf (1235)^ 
6?(1235)tt° 


7.9 
9.2 


± 
± 


1.1 
1.1 


10- 
10- 


-3 
-3 


Bizzarri (1969) a 
Amsler (1993a) a 


ag(1320)w 


1.70 


± 


0.15 


10" 


-2 


Amsler (1994c) 


K*(892)K 


7.05 


± 


0.90 


10- 


-4 


TT ± K^K L 


(%) 


1.05 


± 


0.08 


10- 


-3 


7T°K L K L 




1.5 


± 


0.3 


10- 


-3 


Conforto (1967) 


K*(892)K 


2.11 


± 


0.28 


10" 


-3 






2.70 


± 


0.37 


10- 


-3 


ir°K s K L b 




> 2.51 


± 


0.22 


10- 


-3 


Conforto (1967) 
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Figure 1: Pion multiplicity distribution for pp annihilation at rest in liquid hydrogen. Open 
squares: statistical distribution; full circles: data; open circles: estimates from Guesquiere 
(1974). The curve is a Gauss fit assuming (N) = 5. 



Figure 2: Fraction fp of P-wave annihilation as a function of hydrogen density (curve). The 
dots with error bars give the results from one particular optical model (Dover and Richard, 
1980) using two-body branching ratios (adapted from Batty (1996)). 



Figure 3: The Crystal Barrel detector. 1,2 - yoke, 3 - coil, 4 - CsI(Tl) barrel, 5 - JDC, 6 - 
PWC's, 7 - LH 2 target. 



Figure 4: The silicon vertex detector. 1 - microstrip detectors, 2 - hybrids, 3 - readout elec- 
tronics, 4 - cooling ring (from Regenfus (1997)). 



61 



Figure 5: 27 invariant mass distribution for a sample of 47 events (6 entries/event). 



Figure 6: 7r°7 momentum distribution in pp — > 47r°7 (4 entries/event). The peak is due to 
pp — > ujt\. The inset shows the u;-region and a fit (Gaussian and polynomial background). 



Figure 7: 7r + 7r 7r° invariant mass distribution for tx + ti 7i°i] events. The peak is due to pp — > urj. 
The inset shows the background subtracted angular distribution in the a;- rest frame (see text). 



Figure 8: 7 angular distribution in the cu rest frame for ut](uj — > 7r°7) (see text). 
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Figure 9: Annihilation graph A and rearrangement graph R for pp annihilation into two mesons. 



Figure 10: Following VDM, radiative annihilation can be described by a superposition of two 
isospin amplitudes with unknown relative phase [3. X stands for any neutral meson. 



Figure 11: Energy deposits in the barrel versus polar angle G and azimuthal angle $ for a n 7 
event. The two 7's from 7r° decay cluster near the minimum opening angle (16.5°). 



Figure 12: Energy distribution (24,503 events) of the single 7 in the missing tt° rest frame for 
events satisfying the kinematics pp — ■> 2ir° and a missing it . The full line is a fit. The dotted 
curve shows the expected signal for a branching ratio of 5 x 10~ 4 . 



Figure 13: 90 % confidence level upper limits for radiative pseudoscalar decays as a function of 
missing mass. 
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Figure 14: Squared matrix element for r] — > 3ir°. The straight line shows the fit according to 
Eq. 



Figure 15: (a) u — > 37 Dalitz plot for tt°uj events (62,853 events, 6 entries/event); (b) u — > 37 
Dalitz plot for r/u events (54,865 events, 6 entries/event). 



Figure 16: Dalitz plot of the final state KsKltc (2,834 events). 



Figure 17: Ratio of to c<j production in low energy annihilation. The measured branching 
ratios have been divided by the factor W (Eq. (p0|)). The expectation from the OZI rule using 
the quadratic mass formula (4.2 x 10~ 3 ) is shown by the horizontal line. 



Figure 18: (a,b): OZI allowed production with ss pairs in the nucleon. In (c) an intermediate 
four-quark state is excited below threshold. The production of </> mesons can also be enhanced 
by final state rescattering, K*K — > ir<f) or pp — > 7T0 (d). 
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Figure 19: r/n and KK mass distributions for the ao(980) resonance in pp — > r/nX and KKX 
(in arbitrary units and assuming that no other resonance is produced in these channels). The 
dashed line shows the r/n mass distribution for the same width T' in the absence of KK coupling 
(92 = 0). 



Figure 20: Invariant mass distributions for it !]!]; (a) r/rj mass distribution for 67 events showing 
the two new scalar mesons; (b) 7777 mass distribution for one rj decaying to 3tt° (IO7 final state, 
not corrected for acceptance); (c) 7i r] mass distribution for 67 events (2 entries/event) showing 
the a (980). The solid lines in (a) and (c) represent the best fit described in Amsler (1992c). 



Figure 21: Dalitz plots of 3-pseudoscalar channels. Red and blue regions correspond to high, 
respectively low, event densities; (a) -K r]T] (198,000 events). The Dalitz plot is symmetrized 
across the main diagonal; (b) -k 7c i] (symmetrized, 280'000 events); (c): 3ir° (712,000 events). 
Each event is entered six times for symmetry reasons; (d) tt°KlKl (37,358 events). 



Figure 22: 7r°7r° mass projection mpp — > 3n° (3 entries/ event) with the fit (solid line) described 
in the text. 
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Figure 23: Argand diagram of the tttt scattering amplitude T obtained from a common fit to 
production and scattering data (from Spanier, 1994). 



Figure 24: Isoscalar S-wave production intensities |T| 2 in 3ir° (full curve), 2ir°r] (dashed curve) 
and 2?77r (dotted curve) before multiplying by the phase space factor p. The vertical scale is 
arbitrary (from Spanier, 1994). 



Figure 25: r/r]' mass projection in ir°r)r)'. The full curve is the fit to the 67 data with a scalar 
resonance close to threshold and the dashed curve shows the expected phase space distribution. 
The inset shows the 777/ mass distribution from 7r°r7( — > 2j)t]'(— > r]n + n^). 



Figure 26: 27r°r] mass distribution recoiling against 7r°7r° for events with 10 reconstructed 
photons (6 entries/event). 



Figure 27: rj'ii mass projection in 7r°7r ?/ — > 67 for data (dots) and fits (histograms); (a) (tttt) 
S-wave and a 2 (1320); (b) including a (1450) (best fit). 



66 



Figure 28: The phase shift 5 in elastic Kir scattering from Aston (1988b). The curve shows 
the fit using prescription ( |106| ). 



Figure 29: KK mass projection in ^K^K^. The dashed line shows the fit with one scalar 
resonance, the full line the fit with two scalar resonances. The peak on the left is due to 
/2(1270) and 0,2(1320) and the peak on the right to i^*(892) reflections. The central peak is 
due to interference from various amplitudes. 



Figure 30: Branching ratio for / (1370) and / (1500) decay into KK as a function of a (1450) 
contribution to tt°KlKl (from Dombrowski (1996)). 



Figure 31: dE/dx distribution in the jet drift chamber for 2-prong events with a missing K^. 
The curve shows the expected (Bethe-Bloch) dependence. 



Figure 32: ■k ± K^K l Dalitz plot (11,373 events). 
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Figure 33: x 2 dependence on the fractional contribution from ao(1450) to ^K^Kl- 



Figure 34: Tangent of the nonet mixing angle a as a function of Ri (solid), R 2 (dotted) and 
i?3 (dashed curve). The shaded areas show the experimentally allowed regions for / (1500), 
assuming that this state is qq. 



Figure 35: 7r°7r°r/ (a) and n + n r] (b) mass distributions in pp annihilation at rest into 
7r + 7r~7r°7r°?7, showing the rj and 77(1410) signals. The dashed line shows the result of the 
fit. 



Figure 36: Angular distribution of a (980) ± in the 77(1410) rest frame. The data are shown 
with error bars. The full curve shows the fit for a 0~ + state and the dashed curve the prediction 
for a 1 ++ state produced with the same intensity. 
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